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Abstract 



In the present work we discuss non-linear physics problems such as Nielsen- 
Olesen strings, superconducting bosonic straight strings and static vortex 
rings. We start with a toy model. We search for antiperiodic solitons of the 
Goldstone model on a circle. Such models provide the basis as well as useful 
hints for further research on three-dimensional more realistic problems. We 
proceed with a full research on a U(l) model which admits stable straight 
string solutions in a small, numerically determined area. That model has a 
Ginzburg-Landau potential with a cubic term added to it and can be found 
in condensed matter problems as well. The next part of our research, has 
to do with a 1/(1) x 17(1) model which is the main subject of our interest. 
There, we search for stable axially symmetric solutions which are solitons, 
which can represent particles, the mass of which is of the order of TeV. 
The confirmation or rejection of the existence of those defects is of great 
interest if we consider that LHC will work in the same energy range. In our 
study, we find out that due to current quenching, these vortex rings seem 
to be unstable. We also extend the model of vortex rings by adding higher 
derivative terms which are favorable for stability. After the extensive analysis 
we performed, we conclude that these objects don't seem to be stable. The 
reasons which brought us to this conclusion are explained. 
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Chapter 1 
Introduction 



In order to study topological defects, in general we need to define the order 
parameter, to identify the pattern of broken symmetries and to classify the 
defects we have. A defect is defined as a space region containing the dis- 
continuity of the order parameter. In field theory, the order parameter is 
represented by scalar fields. The topological defect is the singularity in the 
continuum that cannot be removed by continuous deformation of the local 
region in the vicinity of the singularity. Strings, monopoles, textures, domain 
walls or combination of them are some examples of defects. Strings especially 
are of interest to us. 

Current theories of particle physics place chronologically the formation of 
topological defects on the violent stages of the early evolution of the universe. 
If somehow their existence is observationally proved, that would be of great 
importance for particle physics and cosmology, as the notion of topological 
defects exists in most interesting models of high energy physics trying to 
describe the early universe. Those defects would be a direct consequence of 
events of the very early universe. 

The same notion is also present in condensed matter physics. Simple 
examples are the domains in a ferromagnet or the magnetic vortices in thin 
ferromagnetic films [1] or the vortices in Ginzburg-Landau superconductors 
[2]. Domains are regions in which the magnetic dipoles are aligned, separated 
by domain walls. Liquid crystals exhibit an array of topological defects, such 
as strings and monopoles. Defects can also been found in biochemistry, no- 
tably in the process of protein folding. A more specific example is that 
when liquid helium is quickly pressure-quenched around a critical tempera- 
ture, namely T c , Zurek argued that vortices carrying rotation in quantized 
amounts form and represent defects analogous to cosmic strings. More re- 
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Figure 1.1: String intercommutation in a nematic liquid crystal. Further details 
can be found in |llj . 

cently, the existence of vorton solutions is involved in sectors such as QCD [3] 
and high-T c superconductivity [I]. There are differences from the cosmolog- 
ical case where relativistic dynamics must be used and gravity is important, 
but their formation can offer interesting hints for cosmology [3, [6J, [T] . One 
can also think of Bose-Einstein condensates (BEC) (i.e. see [Sj), vortices in 
superfluid Helium-3 and Helium-4 P, [10] , or nematic liquid crystals [HI [12] 
(see also figs JTTlOl) . 

Interest in such defects appears in many branches of theoretical and ex- 
perimental physics such as particle physics models, experimental high energy 
physics, cosmology, observational astrophysics, experimental condensed mat- 
ter physics and, recently, superstring theory. In contrast to other branches 
of physics, it is established that no such topologically stable defects can ex- 
ist in the Standard Model (SM) of Particle Physics. But, as pointed out in 
[13], there is the possibility of quasi-topological metastable defects in popular 
extensions of the SM such as the two-Higgs Standard Model (2HSM) [Ti] . 
This was proved to be the case for codimension one and two defects [H] , but 
the search for particle-like solitons of this kind was not so far succesfull [15] . 
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Such solitons, if they exist, in the electroweak theory are very interesting. 
They will have a mass of a few TeV, and should soon be produced in the 
LHC at CERN. Search for such objects is the main topic of the present work. 



1.1 Formation of defects: The Kibble mech- 
anism 

The explanation of the formation of topological defects, uses a very im- 
portant notion of particle physics which is that of spontaneous symmetry 
breaking. This phenomenon appears when the ground state of the system is 
characterised by a non-zero expectation value of the Higgs field and does not 
exhibit the full symmetry of the Lagrangian. It is believed that the way to 
unification of all fundamental natural forces involves this notion. The Higgs 
field, named after the British physicist Peter Higgs, is a postulated quantum 
field, mediated by the Higgs boson, which is believed to permeate the entire 
universe. Its presence is said to be required in order to explain the large 
difference in mass between those particles which mediate weak interactions 
(the W and Z bosons) and that which mediates electromagnetic interactions 
(the photon). 

This is a spin field which signals the breaking of some symmetry in 
our theory. This means that when our system goes from a state of higher 
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Figure 1.3: Higgs field space and physical space. Above T c the potential has a 
minimum at <j)\ = = (fa, but below T c the new form of the potential has its 
minima on a circle. If m is an element of the coset space G/H and and two 
different states at that circle, then = \&2- Generally, every element of the 

coset space G/H can produce all the other minimum energy states when acting on 
any of them. Thus, G/H = Ai where M stands for the vacuum manifold. On the 
right side of the figure one can have a picture of how a cosmic string forms when 
different regions of space met, having their Higgs field phases as above. Picture 
from |astro-ph/0303504| (A.Gangui). 

symmetry to a final state where it obeys a smaller group of symmetry rules, 
this field represents the symmetry breaking order parameter which acquires 
a non-zero vacuum expectation value. In the stages we described above, this 
field can settle into different ground states or to be more specific, the phase 
of the Higgs field can acquire different values. Topological defects are objects 
that locally restore the original symmetry. More details about spontaneous 
symmetry breaking will be given in the next chapter. 

In the context of the standard Big Bang theory, the spontaneous breaking 
of fundamental symmetries is realized through phase transitions in the early 
universe. The time needed for such a transition is small compared to the 
expansion time. It is believed that all forces were united in the early stages 
of our cosmos and that these phase transitions play a significant role in order 
to end up to what we today know as four different natural forces. It is also 
believed that these transitions are the source of defect formation and if that's 
true, then there is an important link between particle physics and cosmology. 

Kibble in 1976 [16J first noted that, in an expanding universe where are 
separated regions that have no "communication" amongst themselves, due 
to lack of causal contact, we can assume that they have different values of 
the phase of the Higgs field. As the universe expands and cools down, these 
regions expand and some of them come in contact with others. Thus, it is 
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3- DIMENSIONAL 

PHYSICAL SPACE PHYSICAL SPACE 




not impossible that by following a closed path through these regions that are 
in contact, their phases (orientation of the Higgs field) vary from to 2nir, 
where n ^ is an integer and stands for the winding number of the field. 
Then, the "meeting point" between these "bubbles" is the core of a cosmic 
string (figs ll.3l and 11.41) . In other words, a cosmic string forms when, for 
purely topological reasons, the orientation of the arrows cannot be adjusted 
in order to keep the Higgs field in the minimum energy state everywhere. 



1.2 Classification of defects 



The different kinds of defects are due to different symmetries of the Higgs 
field itself. For example, the Higgs field may have two possible states to 
fall into. Then, when regions that have taken opposite choices meet, their 
boundary is a domain wall. This narrow region of the defect is called false 
vacuum and should not exist but it does so, due to the geometry of the Higgs 
field. Now, if Higgs field has circular symmetry it can be represented by a 
"mexican-hat" analogy because ground state can lie anywhere within a circle 
in field space and its position can be denoted by an arrow. This is the case 
shown in the previous section and the defect involved there was a cosmic 
string. 

Consider a temperature-dependent potential which at high temperatures 
has a single minimum at <p c = 0, while under a critical temperature T c it 
transforms into a potential which has a minimum for a non-trivial value of 
the field, <p c ^ (fig ll.5l) . When T < T c the system will try to minimize 
its energy and will move at the new minimum. When this is done, and the 
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T>>T r 



T=0 




Figure 1.5: Temperature dependent potential. Here is an example of such po- 
tential and how it changes shape as temperature falls. It is clearly visible that 
above T c the system has one stable minimum at <p c = 0. As T — > the potential 
develops a minimum at (f) c ^ and when the system settles there, we say that the 
initial symmetry of the system is broken. The ground state of the model does not 



respect the original symmetry of the Lagrangian. Picture from astro-ph/0303504 
(A.Gangui). 



phase transition has finished, our system is ruled by a smaller symmetry 
group, say H. Also name the previous, bigger symmetry group as G. Then 
the coset space G/H produces the vacuum manifold M. and the topology 
of the latter determines the type of defect that will finally arise. Homotopy 
theory helps us by telling us how to map M. into physical space (fig Jl.31) . For 
example 7r\(A4) ^ 1 indicates that the first homotopy group is not trivial 
and practically informs us for the existence of non contractible loops in A4 , 
that is to say cosmic strings. Thus, in models where the symmetry of the 
ground state is associated with a non-trivial homotopy group, topological 
defects exist. 
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1.3 Superconducting cosmic strings 

Cosmic strings are linear vortex defects predicted to have formed at a cosmo- 
logical phase transition during which the vacuum manifold was not simply 
connected. These strings have enormous energy per unit length, namely \i. 
Roughly speaking, for a string created at a phase transition characterized by 
temperature T c , then \i ~ T 2 . For GUT strings, fi ~ 10 22 g/cm. Their grav- 
itational effects though, are negligible. The strength of their gravitational 
interaction is given in terms of the dimensionless quantity Gfi ~ (T c /Mp) 2 , 
with G the Newton's constant and Mp the Planck mass. There are other 
significant astrophysical effects when, for example, strings have supercurrent. 
The latter can make them transform into thin superconducting wires with 
critical current ~ T c thus, interacting strongly with magnetic fields. 

It is useful to note here that the notion of superconductivity is also im- 
portant in our research presented in Part II, where one can find details. 
This phenomenon was first discovered by Onnes in 1911 and exhibits many 
astonishing and interesting features, two very well known of which are: 

• The superconducting material 1 while being below a critical tempera- 
ture T c , allows current flow in its interior with zero resistance. 

• Magnetic fields generally can not penetrate in the interior of a super- 
conducting material, the well known Mei(3ner effect. 2 

There are two kinds of superconductors: 

• Type I: Below a critical field H C (T) there is no penetration of flux 
inside the superconducting material, while above H C (T) the field pen- 
etrates perfectly in the material which has turned into its normal non- 
superconducting state. 

• Type II: Below a critical field H a (T) again there is no penetration. 
The difference here is the existence of a second critical field Ht,(T) 
above which normal state is restored but for values of H C (T) where 

'"Not all elements or combinations of elements can be superconducting. In the years 
passed since its discovery, superconductivity appears in many other materials and the 
challenge for the years to come, is the construction of a material appearing superconducting 
properties below a T c which will be comparable to every day temperatures, say around 
280 °K. We are at - 140 °K for now. 

2 This effect is of crucial importance for the stability of vortex rings which will be studied 
later on. 
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H b (T) > H C {T) > H a (T) is valid, the applied field penetrates par- 
tially in the material which is in a mixed state where both normal and 
superconducting regions coexist. 

1.4 Brief history of research on defects 

• The beginning... 

Skyrme in 1961 [T7], presented the first three dimensional topological defect 
solution arising in a non-linear field theory and proposed that such solutions 
would have an important role in particle physics. Nambu in 1966, suggested 
the cosmological significance of defects. Before such proposals could be taken 
seriously, one had to establish the existence of stable topological defect so- 
lutions in realistic renormalizable theories. An important step towards that 
direction, was the discovery of defect solutions in Higgs and Yang-Mills the- 
ories, notably the Nielsen- Olesen vortex in 1973 [IE], the 't Hooft-Polyakov 
monopole in 1974 [Tj5] as well as GUT monopoles (i.e [2"U]). 

The cosmological implications of symmetry breaking were pointed out by 
Kirzhnits in 1972 [2T], who suggested that spontaneously broken symmetries 
in a field theory can be restored at high temperatures just as they are in 
condensed matter systems. Zel'dovich, Kozbarev and Okun in 1974 [2*2*] . 
argued that domain walls would form at phase transitions in which a discrete 
symmetry was broken. Weinberg in 1974 [23], noted that cosmic domains 
could form and Everett [21], in the same year studied possible interactions 
of domain walls with matter. It was then the pioneering paper of Kibble "Tfj] 
in 1976, which showed that the existence of defects, depends on the topology 
of the vacuum manifold classifying them using homotory theory. 

The circumstances above, led many researchers to start studying those 
defects and their role in cosmology. A new "golden" era for topological 
defects was starting. 

• The "booming"... 

The first interest in studying them comes from the fact that since a typi- 
cal GUT predicts a few phase transitions and because the vacuum structure 
needed to form strings is generically realized, one can assume, also by follow- 
ing Vilenkin, that cosmic strings have a considerable existence probability. 
Cosmic strings appear in many papers especially in the 80's and early 90's as 
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Figure 1.6: Apart from the intersection of two different strings (a), self- 
intersection of a cosmic string can be another mechanism of loop production (b). 
Picture from [135] . 

having important cosmological implications in structure formation by pro- 
viding density fluctuations and thus, becoming the seed of galaxies [T6"|-[3"T]. 
Cosmic strings appear either as infinite straight strings or as rings. How- 
ever, due to disagreement with the observations of the Cosmic Microwave 
Background, this theory was left aside until nowadays, as superstring the- 
ory discusses the possibility for the existence of cosmic superstrings in the 
framework of braneworld cosmologies. Such superstrings could play, under 
conditions, the role of cosmic strings. There was also an observational "ad- 
venture" as we will explain below. 

There are numerous papers on the consequences of a possible existence 
of superconducting vortex rings [32]- [50]. String loops are predicted to have 
formed during phase transitions at the early universe. [15] refers to GUT 
mass scale U(l) symmetry breaking occurence which gives rise to strings 
having important role in many popular cosmological scenarios seen mostly 
in the 80's. The cosmological role of loops has been analyzed in many papers 
such as [51]- [60], where the scenario is that galaxies condensate around oscil- 
lating loops of mass ~ 10 9 M Q . By the time those strings have radiated away 
their energy, the matter density fluctuations are large enough to grow inde- 
pendently. Additional loops can be produced by self-intersection of a string 
or by intersection of different strings (fig ll.6l) . Closed loops are doomed to 
extinction due to oscillation and gradual loss of energy through gravitational 
radiation [25], [26], unless perhaps they reintersect longer strings and become 
reconnected [32]. On the other hand, it's crucial for closed loops to have 
possibility more than 1/2 to intersect themselves and break up into smaller 
pieces in order not to dominate the energy density of a radiation-dominated 
universe [HI]. Qualitative as well as quantitative details were revealed by 
simulations [61]- [67] of a cosmic string network consisting of infinite straight 
strings and a small percentage of loops. 
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• Today... 

Until today, physicists search for a possible signal of a topological defect from 
the early universe through the ways the theories give them for detecting them 
|54j . [68]-[72]. Especially for the case of a ring with supercurrent, when cur- 
rent saturates, such a string will emit particles copiously and may be seen 
as an X-ray or 7-ray source and it is believed to give contributions to X-ray 
background and high-energy (10 20 eV) cosmic rays [71]. In [7_6]-[78] one can 
find observations which were believed to involve, for the first time, a topolog- 
ical defect through gravitational lensing while in [79] the lensing from dark 
matter was excluded. Later observations and analysis though, showed that 
this is not in fact the case [HO]- In fact, it is a very rare case of almost identi- 
cal (up to 99.96% !) elliptical galaxies, separated by a misleadingly small, as 
it proved, distance. In [56] Kibble supports the cosmic string cosmological 
scenario using the above observations and predictions of fundamental string 
theory. Concerning the latter, in [STL 182] one can find details about the pos- 
sibility of the existence of cosmic superstrings, that is to say superstrings of 
cosmic length, and observation of them as cosmic strings. As it concerns the 
observations, if we ever detect a topological defect, that would be a discovery 
of great importance. But even an observed absence of topological defects, 
is very useful too (i.e. [83J), since it imposes constraints on particle physics 
model building. For example, the non-abundance of magnetic monopoles 
inspired the inflationary revolution in cosmology and GUT's models are con- 
strained to provide the requisite amount of inflation. Reviews on all these 
matters can be found in [129j - [TM] and relevant books are [1051 1136] . 

1.5 Solitons 

Topological soliton in general, is a solution of a system of partial differential 
equations or of a quantum field theory that can be proven to exist because 
the boundary conditions entail the existence of homotopically distinct solu- 
tions. Typically, this occurs because the boundary on which the boundary 
conditions are specified, has a non-trivial homotopy group which is preserved 
by differential equations. The solutions to the differential equations are then 
topologically distinct, and are classified by their homotopy class. It is not 
easy to define precisely what a soliton is. Drazin and Johnson (1989) describe 
solitons as solutions of nonlinear differential equations which 

• represent waves of permanent form 
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• are localised, so that they decay or approach a constant at infinity 

• can interact strongly with other solitons, but they emerge from the 
collision unchanged apart from a phase shift. 

There are many equations of mathematical physics which have solutions 
of the soliton type. Correspondingly, the phenomena which they describe, 
be it the motion of waves in shallow water or in an ionized plasma, exhibit 
solitons. The first observation of this kind of wave was made in 1834 by 
John Scott Russell, who followed on horseback a soliton propagating in the 
windings of a channel. In 1895, D. J. Korteweg and H. de Vries proposed an 
equation for the motion of waves in shallow waters which possesses soliton 
solutions, and thus established a mathematical basis for the study of the phe- 
nomenon. Interest in the subject, however, lay dormant for many years, and 
the major body of investigations began only in the 1950s. Researches done 
by analytical methods and by numerical methods made possible with the 
advent of computers gradually led to a complete understanding of solitons. 

Eventually, the fact that solitons exhibit particle-like properties, because 
the energy is at any instant confined to a limited region of space, received 
attention, and solitons were proposed as models for elementary particles. 
However, it is difficult to account for all of the properties of known particles 
in terms of solitons. More recently it has been realized that some of the 
quantum fields which are used to describe particles and their interactions 
also have solutions of the soliton type. The solitons would then appear as 
additional particles, and may have escaped experimental detection because 
their masses are much larger than those of known particles. In this context 
the requirement that solitons emerge unchanged from a collision has been 
found too restrictive, and particle theorists have used the term soliton where 
traditionally the term solitary wave would be used. 

We can classify solitons in three sectors, according to the origin of their 
stability: 

• Topological solitons: They can be found in theories where the vac- 
uum manifold is not simply connected and stability comes as a conse- 
quence of a topological conservation law. 

• Non-topological solitons: These solitons are stabilized by the con- 
finement of a conserved charge as a result of a Noether current arising 
from a symmetry of the theory and not from a topological conservation 
law as above 1851 . 
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• Semi-topological solitons: They are characterized by a winding 
number and are classically stable, being local minima of the energy 
functional, but in contrast to genuine topological defects, they can 
tunnel quantum mechanically and decay to the trivial vacuum. They 
were introduced for the first in [T3l [HI [15], EE] • They are the only kind 
of defects that can exist in the SM of Particle Physics. 
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THEORETICAL FRAMEWORK 



In this part, we present two well known models which provide the basis 
for our research. These are, the Nielsen-Olesen vortex [TS] and the straight 
bosonic superconducting string |103] . These models also provide the theoret- 
ical framework on which we are based. The known solutions of these models 
help us to test our numerical algorithms and to reproduce the solutions. Be- 
fore that, we explain the notion of spontaneous symmetry breaking of a global 
and a local symmetry as well as some other useful notions such as penetration 
depth and Meissner effect in superconductors and phase transitions. 
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Cosmic string field theory 



2.1 Spontaneous symmetry breaking 

The purpose of this section is to introduce the notion of spontaneous symme- 
try breaking (SB) which is very important in our research part. Goldstone 
theorem as well as Higgs mechanism are presented. 

2.1.1 Spontaneous breaking of a global symmetry 

Consider the following Lagrangian density with a complex scalar <fi field, 

£ = A4 W- [7(0, </>*) (2.1) 

with 

U((j>, <P*) = m 2 0> + A(0» 2 (2.2) 

This Lagrangian is invariant under the following U(l) global transformation 

-> e ia (f) (2.3) 

where a a constant. The ground state of the theory is obtained by minimizing 
the potential U 

BJJ 

— = m 2 0* + 2A0*(0*0) (2.4) 

• m 2 > 0: In this case we have a minimum at <fi* = <p = 

• m 2 < 0: In this case, <fi = is a maximum and the minimum is located 
at |0| 2 = -m 2 /2A = k 2 

In the second case, the minima lie along the circle \<p\ = k, which form a set 
of degenerate vacua related to each other by rotation. We can express the 
complex field <fi as follows 

<j>{x) = (r(x) + k)e l9(x) (2.5) 

where r and 9 both have vanishing vacuum expectation values. We substitute 
in L and the potential term there becomes 

U = Ar 4 + 4k\r 3 + 4Xk 2 r 2 - A/c 4 (2.6) 

and we have a mass term for the field r, which means that m 2 = 4Xk 2 , while 
9 is a massless field as there is no mass term for it. Thus, the two massive 
fields (real parts of 0) became one massive and one massless field. In general 
there is the following theorem 

• Goldstone theorem: The spontaneous breaking of a continuous global 
symmetry is always accompanied by the appearence of one or more 
massless scalar spin-0 particles known as the Goldstone bosons. 
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2.1.2 Spontaneous breaking of a gauge symmetry 

Now we consider the Lagrangian density (12.1 1) , which we demand to be in- 
variant under the following U(l) local (gauge) transformation 

-> e ia{x) <f) (2.7) 

This results in the introduction of the electromagnetic field through a covari- 
ant derivative and the Lagrangian density becomes 

£ = + \D^\ 2 - U(cf>, 0*) (2.8) 

where \D^\ 2 = + ieA^d" - ieA»)<jf, F^ = d fJ- A l/ - d v A^ while the 
potential is the same as above. Following the previous steps, in the case we 
consider m 2 < 0, the vacuum is at 



— m 
~2X 



x 1/2 
2 x ' 



(2.9) 



We can set, as above, 



<P( X ) = k + Mx) ^ Mx) (2.10) 

and substitute in the Lagrangian which becomes 

£ = -\f, v F^ + l -{d^{) 2 + \{d^ 2 ) 2 + e 2 k 2 A^ - 2X^1 + ■■■ (2.11) 

where the dots represent cubic and quartic terms. There are two massive 
fields, the </>i and A^ which means that the photon became massive. The 02 
field is massless and can also be eliminated through an appropriate choice of 
the vacuum. The above phenomenon is summarized as follows 

• Higgs mechanism: The spontaneous breaking of a gauge symmetry 
makes the massless field to disappear and we end up with a massive 
scalar field representing the Higgs particle, while the gauge field ac- 
quires mass. 

The Higgs particle is believed to have large mass and this is the reason for 
not being discovered yet. LHC, a proton-proton collider will operate in the 
early 2008. It will work at energies until 7 + 7 TeV, something which means 
that there are possibilities for finding that particle. To have answers in this 
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subject is, maybe, the greatest challenge for the years to come. The following 
table summarizes the above cases. 





Goldstone mode 


Higgs mode 


Before SB 


2 massive scalars 


2 massive scalars+1 photon 


After SB 


1 massive+1 massless scalar 


1 massive scalar+1 massive photon 



Finally, if the ground state of the system, obeys the initial symmetry of the 
Lagrangian, the system is said to be in the Wigner mode. 



2.2 Effective potential 

The discussion of symmetry breaking is somewhat simplistic because we use 
classical potentials to determine the expectation value of the Higgs field </> in 
the models above. We have to keep in mind that, in reality is a quantum 
field interacting with itself and with other quantum fields and the classical 
potential V (</>), is modified by radiative corrections. The corrected potential 
for 0, is called effective potential V e ff and can be evaluated perturbatively 
as an expansion in powers of coupling constants 

Veff{<t>) = V{4>) + + V 2 {<f>) + ■■■ (2.12) 

where V(4>) is the classical potential and V n (4>) is the contribution of Feynman 
diagrams with n closed loops. In some models, the radiative corrections are 
negligible, while in others they can alter the character of symmetry breaking. 
One example of the latter is the Higgs model 



£ = -- V" + D^Dy - V{\(j>\) (2.13) 

with V(|0|) = /^ol^l 2 ; where /io a constant. The one- loop contribution to 
V eff (4>) is given by [M] : 

WH^I^TH (2-14) 
where a the renormalization scale. If we introduce a dimensionless quantity 

then, for v > 0.45 the V e //(0) has a single minimum at = 0, but for 
v < 0.45 there is another minimum at a non-zero value of while for v < 
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Figure 2.1: An example of a first order phase transition. The plots are for (3 = 

2 < Pcrit, P = Peru, P = 2.2 > [) crit . Parameter is a = 1. 

0.37 the latter minimum goes deeper than the one at = 0. In the last 
case, the absolute minimum is at a non-zero value of <fi and the symmetry is 
spontaneously broken. 

In Chapter 4, we examine a modified Ginzburg-Landau potential. The 
modification has to do with the addition of a cubic term. Such term results 
from the 1-loop radiative corrections on the quartic potential [102] . 

2.3 Phase transitions 

2.3.1 First order phase transition 

In order to observe phase transitions and to find the distinguishing features 
we suppose appropriate effective potentials. For a first order phase transition, 
consider the following potential which we examine in detail in Chapter 4: 

^(I0l) = |l0| 2 ^l0| 2 -fl0l + ^ (2-16) 

This potential has an obvious minimum at |0| = 0. This corresponds to the 
symmetric phase of the system. For f3 > 2, it develops a secondary minimum 
at |0| = 1. The critical value for this potential is (3 = 3/\/2 = Peru- This 
happens because for P > P cr it, the secondary minimum goes deeper than the 
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((0.5-6)*x*x/1 2)+(0.25"x*x*x , x) 
(0 25*x*x*x*x) 
((12-6)*x*x/1 2)+(o'25'x*x*x*x) 



Figure 2.2: An example of a second order phase transition. The plots are for 
T = \/0.5 < T cr i tl T = T cr u, T = \/l2 > T cr i t . Parameters are rj = 1, A = 1. 

one at \(f)\ =0. When (3 = f3 cr i t , the secondary minimum makes the potential 
zero, just like the original one at \<fr\ = (fig J2.ll) . The non-zero expectation 
value of \d>\ is 



0+ VP 2 



(2.17) 



Generally, in a first-order phase transition, a new broken symmetry phase 
usually nucleates after the temperature falls some degrees below the critical 
temperature T c . Separated areas of the new phase form independently and 
expand, resulting in a local selection of the broken symmetry vacuum. As 
a result of these independent selections, the final configuration may not be 
able to get rid of the locked-out fragments of the original vacuum (before 
transition). 



2.3.2 Second order phase transition 

In order to observe a second order phase transition, consider the following 
temperature-dependent potential: 



with 



f/(|0|,T) = m 2 (T)|0| 2 + ^| 4 



™ 2 (T) = A(T 2 -6rr 



(2.18) 
(2.19) 
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The critical value for this potential is T crit = ?7v^6- Above this temperature, 
the mass-squared term becomes positive, thus the potential is a sum of posi- 
tive terms with a unique minimum at \(f)\ = which represents the symmetric 
phase (fig i2.2l) . When T < T cr i t , then m 2 (T) < and the order parameter 
1 4>\ acquires a non-zero expectation value which is 

\ ( t ) \ = ^{ T cru-T) 1/2 (2.20) 

In second-order phase transitions, the phase transformation occurs almost 
simultaneously throughout the volume. However, unless the critical temper- 
ature is traversed infinitesimally slowly, the resulting broken symmetry phase 
will contain many distinct regions with different choices of the local vacuum. 
This is because the selection of the new vacuum has to be, somehow, com- 
municated if the same choice is to made elsewhere, and the speed at which 
this transformation can be propagated is finite. 

The main difference between these phase transitions is that, at the first 
order phase transition the symmetric phase (|0| = 0) remains metastable 
when passing the critical value f3 cr u and is often called "false vacuum". A 
defining feature of the second order phase transition, is that the order pa- 
rameter |0| grows continuously from zero, as the temperature is decreased 
below T crit . 

2.4 Superconductivity and penetration depth 

Ginzburg-Landau theory is used to model superconductivity. It examines 
the macroscopic properties of a superconductor with the aid of general ther- 
modynamic arguments. Based on Landau's previously-established theory of 
second-order phase transitions, Landau and Ginzburg argued that the free 
energy F of a superconductor near the superconducting transition can be 
expressed in terms of a complex order parameter which describes how 
deep into the superconducting phase the system is. The free energy has the 
form 

F = F n + a\^\ 2 + ||*| 4 + -L|(-iW - 2eA)*| 2 + ( 2 .21) 
2 2m 2/io 

where F n is the free energy in the normal phase, a and (5 are phenomenological 
parameters, A is the electromagnetic vector potential, and B is the magnetic 
field. By minimizing the free energy with respect to fluctuations in the or- 
der parameter and the vector potential, one arrives at the Ginzburg-Landau 
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equations. These equations produce many interesting and valid results. Per- 
haps the most important of these is its prediction of the existence of two 
characteristic lengths in a superconductor. The first is a coherence length £, 
given by 

f=A/-^T-r (2-22) 
^ y 2m\a\ V 7 

which describes the size of thermodynamic fluctuations in the superconduct- 
ing phase. The second is the penetration depth A, given by 



where is the equilibrium value of the order parameter in the absence of an 
electromagnetic field. The penetration depth describes the length to which 
an external magnetic field can penetrate the superconductor. In general, if 
one considers a superconducting semi-space at x > 0, and weak external 
magnetic field Bq applied along ^-direction in the empty space x < 0, then 
inside the superconductor the magnetic field is given by 

B(x) = B e~ x/X (2.24) 

In atomic units, one can see that the penetration depth is, in fact, 

A oc — — (2.25) 
e^ 

We will use this result in Chapter 5. 

The formula f)2.24p . can be derived as follows. Suppose that an electric 
field E momentarily arises within a superconductor. The superconducting 
electrons will be freely accelerated without dissipation and their mean veloc- 
ity will be 

m— = -eE (2.26) 
dt v 

The current density is 

j = -env (2.27) 

where n the number of electrons per unit volume, and the above formula, 
with the help of (12.261) can be written as 

f = — E (2.28) 
dt m v ' 
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x = 



Figure 2.3: The magnetic field drops exponentially inside the superconductor. 



If one substitutes into Faraday's law, this gives 



„ „ ldB „ . 

V x E = V x j 

c at 



ne 



B 



mc 



Plug the result into Maxwell's equation and this results to 



47r. _o„ 47me 2 „ 
13 - — B 



V x B = — j V'B 



Now, rename mc 2 /4irne 2 = X 2 and you get 



(2.29) 



(2.30) 



V 2 B = -B 



(2.31) 



where A is the London penetration depth. The last equation can be easily 
solved to give the result in (12. 24ft . 
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2.5 Nielsen-Olesen vortex 
2.5.1 Introduction 

The simplest theory exhibiting string solutions is that of a complex scalar 
field 4>(x), in a model described by the following Lagrangian density: 

C = d^*d^-V^), V( ( j ) ) = ^\ 2 -V 2 ) 2 (2.32) 

which has a global U(l) symmetry under the transformation <ft — > <pe ia where 
a a constant. The condensed matter analogue of global strings is the vortices 
in superfluid Helium-4 [86] . 

If we move a bid further and consider a gauge symmetry under the trans- 
formation — > (f>e m ^ 1 ' then we are led to the introduction of a vector field 
Afj, while the replacement of by the covariant derivative = <9 M — 2eA M is 
required. Due to the existence of the vector field A^, a relevant "magnetic" 
field exists as well. This theory describes the Nielsen-Olesen string [TH] . 



2.5.2 The model 

Consider an Abelian Higgs model in two dimensions with a charged complex 
scalar field (f> and gauge field A 11 . The Lagrangian density of the model 
under consideration is: 

C = D^*Dy - \f, u F^ - ^(|0| 2 - rff (2.33) 

where _D M = d^ — ieA^, the covariant derivative and = d^A u — d u A^. The 
energy of a static vortex configuration in two dimensions follows: 

E = jd 2 p (W| 2 + \{E 2 + B 2 ) + V{\<l>\) \ (2.34) 

The vortex energy must be finite which means that field configuration must 
have a form that decreases to fast enough, as p — > oo. This requirement 
imposes some conditions on the asymptotic field configuration which are: 



\<f>\->V , p^oo 
1 d$ 
ep dip 



A v — > — — H , p — > oo 
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where & is the phase of the Higgs field (ft 1 . The last condition on the gauge 
field, can be understood through the requirement that the covariant deriva- 
tive must be zero as p — > oo and 0^0 there. When following a closed 
path around the vortex, the Higgs field must return to its original value 
0(0) = 0(27r) thus 1 = exp 2iirn, which means that the winding number 
n must be an integer. Also, one expects somewhere in the vortex core to 
have 0(x) = for a non-trivial winding number. The fact that the winding 
number is an integer results to quantization of the magnetic flux. 

1 f 2w d$ , e / e f „ 27m , 

n = — —dip = — 6 Adl = — B-dS^$B = 2.35 

2nJ Q dip * 2tt/ s i 2nJ e 1 ; 

where S 1 is a circle of infinite radius centered on this string. This flux 
quantization is a result of the vanishing of the covariant derivative which 
determines A in terms of <9j0. The phase of must change by 27m which 
forces the flux to be quantized. A rough estimation of the size and mass 
of our vortex can be done. These two characteristic features of the vortex 
are determined by the regions over which the scalar and gauge fields start 
to differ from their asymptotic values. These distance scales p s ,p v (s for 
scalar, v for vector), are related to the Higgs and vector particle masses off 
the string. 

p s « m- 1 = (^Xe)' 1 (2.36) 
p, « m~ l = (v^eri)- 1 (2.37) 

Written in terms of these lengths, and assuming p v > p s (type II supercon- 
ductor 2 ), the energy reads: 

E*2^[ln^) + -±- + \r ] *pl] (2.38) 

Ps e 7 Pv J 

where the first term represents the gradient energy of the scalar field, the 
second is due to the fact that the magnetic flux does not wish to be confined 
and the last term is the cost of the difference of |0| from its vacuum value 77. 
Finally, this energy is minimized by replacing the characteristic lengths with 
their relative masses as above. Then we are led to the mass of the vortex 
configuration: 

/i^27mW— , (2.39) 



1 Suppose 4>{p) = f(p)e me , 9 the phase of the Higgs field, n the winding number. 

2 Our model can be compared with that of Ginzburg-Landau which has the same La- 
grangian with ours and p s can be thought as the coherence length -having to do with 
the transition layer from superconducting to normal state- and p v the London penetra- 
tion depth -referring to the exponential decay of the magnetic field at the surface of the 
superconductor- in terms of superconductivity. 
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a formula that in three dimensions gives the mass per unit length of the 
string. 

The vortices in the Abelian Higgs model have condensed matter ana- 
logues like flux tubes in superconductors |124j but there are differences |125l 
126J because Nielsen- Olesen vortices exist in a vacuum background while 
superconductor vortices are amongst charged bosons called Cooper pairs. 



2.5.3 Ansatz and equations 

The Euler-Lagrange equations of our model are: 

{d^-ieA^){d» -ieA^+^iW-rf) = (2.40) 

d,F^ = f (2.41) 

where j u = 2elm[<j)(d u —ieA u )(f)]. We work in the Lorentz gauge (<9 M A M = 0) 
and we rescaled the fields (0 — > ?7 _1 , — > rj^^^A^ and x — > rjx), while we 
also set r) = 1 for convenience. A proper ansatz for a string lying on the 
z-axis would be: 

m = f(py nv (2.42) 

71 

A(p) = -a{p)<p (2.43) 
ep 

where p, (p are the cylindrical unit vectors. We use cylindrical coordinates 
(p,<p,z). According to the above, the appropriate asymptotic conditions 
while p — > oo and boundary conditions at p — > are: 

/(p^oo)^l , a(p ^oo)^l (2.44) 
/(p = 0)=0 , a(p = 0)=0 (2.45) 

With the above ansatz, we are led to the following system of non-linear 
differential equations: 

d 2 f ldf n 2 f 2 A 2 

d^ + -pdp-^ {a - 1} "2 /(/ " 1)=0 (2 ' 46) 

^-^-2eV 2 (a-D = (2.47) 
dp z p dp 

It is helpful to acquire an approximate asymptotic solution as p — ► oo. This 
is: 

a(p) w 1 - C( v ^exp(-v / 2ep)) (2.48) 
/(p) w l-(9(exp(-v / Ap)) (2.49) 
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Finally, another way to acquire the solution of the system is through the 
minimization of the energy functional which is 



E = 2ti 



palp 



(d p a) 
4p 2 



[2.50) 



2.5.4 Virial theorem 

Before the numerical results, it is important to have a way to check whether 
a final configuration of an algorithm can be accepted or not. 

The energy density for a static solution, which is of interest here, reads: 



-.00 



F^ + |A0| 2 + ^(|0|)=£ 



(2.51) 



The solution of our physical problem, has to minimize the energy of the 
system: 



SE SE 8E 



Define 



SE SE 







SE 

SA, 



■diAj 



which together with 



SE _ dT 00 dT 00 
5$ ~ ~cW ~ k d(d k 1>) 
can be written in a shorter form fi = djGij where 



ds 



(:,J :6 " ' <)(i) j0 ) 



di<, 



de 



-di 



-d,A 



d{d^) d(djA k ) 
This means that any static solution of the field equations satisfies 





djGij 



But, we have 



de 



d(d j( f>) 

de 
d(dj(f)*) 

de 
d(d 3 A k ) 



(DAT 
DA 

djA k — d k Aj = Fj k 
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which leads to 

G .. = eSy-iDjWdifP-DjtpdiP-FjkdiAk 

= £5 i:i - (Dj(j))*Di(f) - Dj(f)(Di(j))* + ieA4{D^)* 

- ieAi(f>* Dj4> - F jk F ik - F jk d k Ai 

= eSij - (Di(j))*Dj(f) - (DrfYDrf - F lk F jk - d k {F jk Ai) 

+ A, (d k F jk - ie{<f>*D j( f> - (Z^) V)) 

where on the second line we just add and subtract the necessary terms in 
order to derive the first four terms of the third line. Now, on the third line 
we have 

djd k (djA k - d k A )A t = (d]d k A k - dldjA^Ai = (2.57) 
and from the field equations, the last term of the third line is also zero as 
5 E 5T 00 

5A =0= 8{d^) = - ie ^* D rf ~ ^ D ^T) + d * F ik = (2.58) 

thus 

Gij = eS i:j - F ik F, jk - (Di<f>)*Dj<f> - {DrfYDrf (2.59) 

but 

Fi k Fj k = e nik B n e m j k B m = B Sy — BjBj (2.60) 

which means that 

= ^ - ^ + \D<j)\ 2 + v\ Sij + B t B 3 - {D i( f>)*D j( j> - (Dj^)*D i( f> (2.61) 

From (I2.56P one obtains 

J d 3 xG ik = = J d 3 xdj(x k Gij) = J dSjXkGij (2.62) 

The last equality follows from Gauss theorem, while the first is due to the 
fact that 

dx k 

dj^XkGij) = -Q^Gij + x k (djGij) = 5 kj Gij = G ik (2.63) 
where we also used (12. 56ft . 
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• Virial relation for Nielsen-Olesen string model 



The formulas above are general and can be applied to any model which is 
being examined for static solutions. For the model at hand, take the trace 
of equation (12.621) 

J d 3 xTrG = J dS j x i G lj = (2.64) 
while equation (I2.6ip leads to 

TrG = -^B 2 + \D(f)\ 2 + 3V (2.65) 

Thus, by taking a cylindrical surface of integration around the infinite string 
with cylindrical symmetry we are led to a virial relation 

d 2 x(^B 2 -V^j = (2.66) 
Using V x A = B and remembering that in our case 

A Vm = \{\^-rnlf (2.67) 
we acquire the formula for the magnetic field 

B = 1^ (2.68) 
ep dp 

and rewrite in a more convenient way for our numerical purposes, the virial 
relation 

1 I da\ A 



2 *'*ISV^j -i(f--™ir-)=° (2.69) 

Numerically though, we "break" the integral into two parts. Consider for 
example 



_ H/1H/2II 

want this index as small as possible. 



Then we must have I\ + I2 = 0. We define the index V = f | + f ^ and we 
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Figure 2.4: Nielsen-Olesen string. The graph on top is the output of a relaxation 
program. The graph in the middle is the result of a Newton-Raphson algorithm 
and the bottom graph is done through energy minimization. All graphs show f(p) 
and a(p), for n = 1, A = 2e 2 = 4. 

2.5.5 Numerical solution 

There are many different numerical methods for solving the above system of 
equations. In figure 12.41 one can see the result of three different algorithms. 
Details concerning the Relaxation algorithm can be found on page 762 of 
[104] and a brief outline in chapter 7 (Part III). Details about the Newton- 
Raphson algorithm can also be found in chapter 7. Finally, as it concerns 
minimization algorithm, many details can be found in chapter 7, while we 
use minimization subroutines from page 425 of [104] . 
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2.6 Straight bosonic superconducting string 
2.6.1 Introduction 

Superconductivity can be understood as a spontaneously broken electromag- 
netic gauge invariance. When the gauge invariance is broken, then any mag- 
netic field applied at the boundary of a superconductor decays exponentially 
towards its interior. This happens due to the screening of the superconduct- 
ing current which flows along the boundary. This is known as Mei/3ner effect. 
Having this mechanism in mind, one can support the idea that cosmic strings 
can become superconductors and behave as thin superconducting wires with 
an enormous critical current. For this to happen, a charged scalar field must 
acquire a non-zero expectation value in the neighborhood of the string core. 
Practically, this points out the need of two complex scalar fields <j) and a in- 
teracting with seperate and U a (i) gauge fields _R M and respectively. 
The breaking of 11^(1) is responsible for the existence of vortices, while if 
we identify U a (l) with electromagnetism, then its breaking inside the string 
will result to supercurrents flowing along the string. If more fields are added 
to the Abelian Higgs model describing an Nielsen-Olesen vortex, then the 
appearence of currents in the core of the defect is possible. The notions in- 
volved in this case, which is one step further from the previous one, will be 
the basis for the case of vortex rings later on. 



2.6.2 The model 

The above description takes us to the following Lagrangian density for 
this model: 

C = {D^l 2 + |ZVI 2 - \B,uB^ - \f, u F^ - \a\) (2.71) 

with potential of the following form: 

ni0l,H) = ^(|0| 2 -^) 2 + ^^| 2 -^ +|^IVI 2 -^ ( 2 - 72 ) 

and D^a = (d^ + ieAJa, = (d^ + igR^fa F^ = d ii A v - d u A^ = 
d^R u — d u R^. 

Away from the string, we have the vacuum of the theory where U^l) is 
broken. There, we have < >= 77 and the potential has the form: 

AJcH 



V 



^H) = ^(H 2 + 2K-m 2 )) (2.73) 
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which takes us to the condition: 

vrf - m 2 > => vrf > m 2 (2.74) 

The latter ensures that away from the string: < o >= 0. 

A topological defect, here a string, is a region where the original bigger 
symmetry of our system is restored, as mentioned in the introduction. There, 
we have < <fi >= and the potential becomes: 

V(\M=0,W\) = ^V 4 + ^W\ 4 -^W\ 2 (2.75) 
Breaking of U a (l) inside the string translates to 3 < a >= -7= and the 

V Act 

requirement of a non-zero potential means another condition: 

Y"' > t (276) 

The Euler-Lagrange equations follow: 

D.D^ + ^(|0| 2 - (r? - fw?)y = (2.77) 

d^o + y (i-r + ^ (i^r - v))" = (2 - 78) 

d^ v = J u = -ig (^D u (f) - 0£>V) (2.79) 
d^ v = J V = -le (oD v o - uD V u^ (2.80) 

Exploiting the equations with the currents above, useful and general rela- 
tions can be extracted which can be helpful in order to write the differential 
equations for the gauge fields and R^. Maxwell equation V x B = J will 
be also used. Suppose a very general ansatz in three dimensions for and a 
such as 

= M x i V, z)exp(iip 4> (x, y, z)) (2.81) 
cr = a Q (x, y, z)exp(iip a (x, y, z)) (2.82) 



3 < a >t^ inside the string and dV ^ d J ) ' <y " 1 = — > a 
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Then we get: 

'd 2 A z 3 2 A Z \ 3 2 A X 3 2 A y 2 . „ . , 

^ + ^) + H + m = J ' = - 2e<T ° M - " a ^ (2 - 84) 

a 2 ^ d 2 AA <9 2 A V d 2 A z 2/ A n , , 

TfcT + -&r j + + = J * = - 2e ^ eA * - d ^ < 2 ' 85) 

# + w) + SI + & = = " < 2 - 86 ' 

a/ + j + + = J > = - 2 ^' 9 ' 4 » " (2 - 87) 

IF + + S + St = * = " 2 ^(^ - < 2 ^> 

The first three equations lead to a differential equation for and the last 
three for R^. 



2.6.3 Ansatz and equations 

The ansatz we use for the fields in the case of an infinite straight supercon- 
ducting string lying along the z-axis, is the following: 

a = a(p) (2.89) 

A, = ^d^z (2.90) 

4 = f{p)*r* (2.9i) 

Rfi = — (2.92) 

where is the polar angle and n the winding number of the field <j). We use 

cylindrical coordinates (p, <p, z). These gauge fields have "magnetic" fields of 
the following form: 

B = -1^0 (2.93) 
e dp 

B R = ±» (2.94) 
99 dp 
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where p, ip, z are the cylindrical unit vectors. In order to have finite energy, 
the function P(p) must vanish for large p while the regularity condition at 
small p demands that .P(O) = 1 and 0(0) = 0. The scalar field a does not 
have dependence on p at the origin or else we would have a singularity there. 
The vector potential describes the effects of any current present in the 
string. For large p the function I(p) oc I to tln(p), where I to t is the total current 
inside the string. 

The natural unit of length is the radius of the string 5 = l/^-y^A^ so we 
choose to rescale our fields: a = Yrj, = Xr], I = Qr) and p — > x = ^JX^qp. 
Thus, the system of field equations to be solved becomes: 

X " + -X' - (P ~ 1)2 X - -X 3 + -X- vXY 2 = (2.95) 
x x 2 2 v ; 

Q" + -Q' - e 2 Y 2 Q = (2.96) 

x 

Y" + -Y' - YQ 2 - XY 3 + rh 2 Y - vX 2 Y = (2.97) 
x 

P" --P' -g 2 X 2 (P -1) =0 (2.98) 

x 

where prime denotes differentiation with respect to x and A = A^/A^,, v = 
v/Xfa g 2 = g 2 /X l j ) , e 2 = e 2 /X l j ) , rh 2 = m 2 /X ( f ) rj 2 . According to these changes 
the conditions (12.741) . (I2.76P become: 

v > rh 2 , A > 2m 4 (2.99) 

Notice that if we approximate X(x) by tanh(x/2) for an ordinary cosmic 
string and search for a solution with vanishing electromagnetic current (Q = 
0) then the F-equation decouples from the others and becomes: 

Y" + -Y' - XY 2 - v tanh 2 (x/2)F + rh 2 Y = (2.100) 
x 

Finally, the energy functional which one has to minimize in order to 
get the final field configuration of this model, is: 



E = 2ti pdp 



oo 







(d p X) 2 + (d p Y) 2 + ^- 2 {d p P) 2 + ^- 2 (d p Q) 2 + 



2r]e 2y " ' 2g 2 p 2 



2 



+ -AQ- ifX 2 + -Y 2 P 2 + ±(X 2 - l) 2 + V 4 + t;Y 2 \Q 2 - — 

p 2K rj 8 V 7 4 2 \ v 
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2.6.4 Numerical solution 

We used Newton- Raphson algorithm (chapter 7) to solve the system of equa- 
tions above with an initial guess that satisfies the boundary conditions and 
asymptotics of the fields. The virial relation for the problem at hand, changes 
slightly because of the existence of the scalar condensate cr(x) and another 
"magnetic" term is added. Thus the virial relation we check, is: 



We also found a solution to this model through minimization algorithm. 
The output of both algorithms is exhibited in figure 12.51 and the small differ- 
ences one may observe are due to small changes in the parameters. A very 
helpful review and analysis concerning especially the Y field, which is the 
charge condensate, can be found in [87] . 




(2.101) 
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Figure 2.5: Straight bosonic superconducting string. Both graphs on top are 
the output of energy minimization algorithm. Parameters are (A, fh 2 , g 2 , v, e) 
= (1,0.65,0.25,0.75,1). Bottom graphs are the output of Newton-Raphson al- 
gorithm. Parameters are (A, fh 2 , g 2 , v, e) =(1,0.62,0.25,0.75,1). Letter a denotes 
the scalar field, b the charge condensate and c the magnetic field associated with 
the a field. 
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PART II 



SEARCH FOR STRINGS AND SOLITONS 

This part includes the research. Briefly, the first section involves a search 
on antiperiodic solitons on S 1 which is a based on [HE]- The second section 
does a full numerical analysis of a straight string with a Ginzburg-Landau 
potential with a cubic term added to it. Such a potential is used in condensed 
matter physics as well [101] . The third section presents a numerical search on 
static bosonic superconducting vortex rings which is based on [r5].[T5] and 
can be seen as a continuation of [JS], where spherically symmetric solutions 
were examined. The final section is based on the model of static vortex rings 
which is now extended by adding higher derivative terms. We analyze the 
behavior of the solutions of the model under these modifications. 
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Chapter 3 

Antiperiodic solitons of the 
Goldstone model on S 1 
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3.1 Introduction 

The purpose of this section is to present the complete list of static classical 
solutions of the Goldstone model on a circle S 1 of radius L in 1 + 1 dimensions 
and the corresponding bifurcation tree together with a study of the stability 
of our solutions. It comes as a supplement to a previous note [HE] which was 
searching for such solutions but with periodicity condition imposed there. 
Jacobi elliptic and standard trigonometric functions are used to express the 
solutions found and stability analysis of the latter is what follows. Classically 
stable quasi-topological solitons are identified. Many of the results have 
obvious similarities with those in [88]. We notice those and make comparisons 
with this reference. We also check the case of mixed boundary conditions and 
place a note with our conclusions in the end. The following analysis has both 
mathematical |108j and physical interest as it can be useful for the search 
of stable solitons in the two-Higgs standard model (2HSM) or the minimal 
supersymmetric standard model. 

3.2 The classical solutions 

The Langrangian density of our model is 

C = \(d,<p 1 ) 2 + ^d,<p 2 ) 2 -V(<j )l ,<j )2 ), fi = 0,l (3.1) 
with (pi and 2 real Higgs fields and potential of the following form: 

V{<k,<h) = \{<p\ + <Pl-l) 2 (3.2) 

We impose antiperiodic boundary conditions on the scalar fields with coor- 
dinates x e [0, 2ttL]. 

</>i (x + 2ttL) = -<f>i (x) , 2 = 1,2 (3.3) 

The general static solution of model (13. ip can be expressed in terms of Jacobi 
elliptic functions: 

0i = VR cosf2, 2 = VR sinQ (3.4) 
R(x) — a± + ci2sn 2 (V2 A(x — Xq) , k) (3.5) 
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n(x) = c l'w) dy 



(3.6) 



where C (a x (k, A) , a 2 (k, A)) ,x ,£ are constants while sn denotes the Jacobi 
elliptic function; sn(z,k), sn 2 (z,k) are periodic functions on the real axis 
with periods 4K(k) and 2K(k), respectively. K(k) is the complete elliptic in- 
tegral of the first kind. Inserting fl3.4l) - fl3.6p into the equations corresponding 
to (13. ip leads to several conditions on the parameters and as a consequence 
to three different types of non-trivial solutions. We start with the simplest, 
the trivial solution. 



3.3 The trivial solution 



The energy functional for static configurations is given by 



E 



2ttL 



dx 



and the field equations are 

d 2 6 



dx 2 



d 2 c 



dx 2 



2 \ dx 



M0? + 02 



I 01 + 02 



2 \ dx 



- 1 

- 1 








V{<f> u <h 



(3.7) 



(3.8) 
(3.9) 



Apart from the vacuum solutions which have E vac = 0, one can immedi- 
ately think of the simplest solution, the trivial one 



, En 



Ltt 



(3.10) 



which exists for all values of L. The corresponding small oscillation eigen- 
modes, labeled by j, have 



u 2 (j) = -((j + 1/2) 2 -L 2 ), j = 0,1,2... 



(3.11) 



This solution is stable until L = 1/2 because o) 2 (0) < for L > 1/2. Many 
additional solutions, some of which were discussed in [13] bifurcate from the 
solution <pi = <p2 = at critical values of L. 
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3.4 Three types of non-trivial solutions 

Now, we present without many details the three non-trivial solutions and 
mention their similarities with those in [88J. The simplest case one can 
think, is to set 02 = so that from (13. 5ft we obtain R(x) = a\ and the 
solution becomes 



I? C ° S L 1 



-^sin(^) (3.13) 

where N is an integer. This solution is called type-I. One can observe that 
the above solution reduces to ( 13 . 1 Of) when L = N + 1/2 i.e. when one of the 
u) of equation (13. lip crosses zero. The Higgs field winds (2N + l)/2 times 
around the top of the Mexican hat. It's energy is given by 

E^N) = 1{n + 1 ^-^{n + \)\ (3.14) 

If we denote the above result which corresponds to the antisymmetric case we 
study here with Ef and the result of [BE] for the same type of solutions with 
Ef where S means "Symmetric" and A "Antisymmetric", then we notice 
that 

E?{N,L) = E?(n+ 1 - j 

Another choice is to set a\ = on ( 13. 5ft so the solution now involves the 
Jacobi elliptic functions. C becomes zero as well and the solution is 

0! = 2kAsn (V2Ax,fc) , 2 = , A 2 = 2 ^ +k2) (3.15) 

This solution is called type-II. In fact, it corresponds to an oscillation of the 
Higgs field in the 02 = plane about the origin <pi = = 2 . If we take 
account of the antiperiodicity condition (13.31) . then the argument k of the 
Jacobi elliptic function is related to the radius L of S l through the following 
formula 

71 
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for some integer m. When we reach the limit of k — > (i.e. L — > m + |) the 
solution (13.1 5p approaches (13.101) . The energy of the solution above is given 
by the integral 



E 



n 



4(2m + 1) 



K(k) 



dy 



k sn (y, k) 



1 + k 



2\ 2 



V2A (1 + k 2 ) 2 
and by means of Elliptic integrals it becomes 



+ 



2k 2 



(3.17) 



E 



n 



4 (2m + 1) 
V2A (1 + k 2 f 



K{k) + 2k 2 - 5) + ^ (k 2 + 1) 



24 v ' 3 

Two specific values of k with the corresponding results follow for k = 

(2m + 1) 7T 



(3.18) 



E n (L = m + 1/2, m) 



(3.19) 



and for k — 1 



(L = oo, m) 



4(2m + i; 
3^2 



(3.20) 



The comparison of Efj which represents our solutions, with E n which repre- 
sents the symmetric case of the same type of solutions studied in [HE] , implies 
that 

1- 



E n [m 



(3.21) 



If none of ai, a<i are zero then we are led to type-Ill solutions where we have 
the following conditions for a 1; 02 and C 2 



ai = ^(l-2A 2 (l + A; 2 )), a 2 = 4k 2 A 2 



(3.22) 



C 2 = 4(l+(4A; 2 -2)A 2 )(l + (4-2A; 2 )A 2 )(l-2A 2 (l + A; 2 )) = 

= ^(l + (4A; 2 -2)A 2 )(l + (4-2A; 2 )A 2 )a 1 (3.23) 
9 

Here, one can explicitly observe the fact that when ai = then C 2 = 
as well. In addition, R and C 2 must not be negative. This means another 
condition 
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A ' 2 £ IfTTP) < 3 - 24 > 

with the equality leading to type-II solutions. 

In order for Q to satisfy Q (x + 2nL) = Q (x) + it and R to be periodic 
(so for the solution to be antiperiodic as we want) on [0, 2nL] the following 
equations must hold: 

Cj j^dy={2n + l)n (3.25) 



o 



(2m + 1) K(k) 

1 (3.26) 



v^vrA 

where m, n are positive integers which respectively determine the number of 
oscillations of the modulus of the Higgs field and the number of the Higgs 
field windings around the origin 0i = = 02 in a period 2txL. 

Now we return to (13.221) and (13.231) and analyze further type-Ill solutions. 
Solving these equations for k = 0, remembering that K(0) = tt/2, we find 
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the critical values of L where these solutions start to exist: 

A 2 = — L 2 = - (2n + l) 2 - — (3 27) 

6 (2n + 1) 2 - 4m 2 4 v ' 2 K ' ' 

The expression for A 2 together with (13.241) leads to the condition 2n + 1 > m 
on the integers m and n (m and n / 0). For any n and m ^ 2n + 1 the 
coefficient of sn 2 in the integral (I3.6P is proportional to k 2 and if one expands 
the solution in powers of k 2 is led to the following formulae 

A 2 = — 2 ( l + ¥)+0 (k*) (3.28) 

6(2n+l) 2 -4m 2 V 2J V ' 

L 2 = 3 -(2n + l) 2 -^ + (k*) (3.29) 



„ 4((2n + l) 2 -m 2 ) 2 (2n + l) 2 

C 2 = -^ J — ^— + (k 4 ). 3.30 

3 ( 3 (2n + l) 2 -2m 2 ) 3 V * 

This expansion is necessary as there is no closed form for integral (13.251) . In 
the limit k = we observe that the 2n solutions of type-Ill yield the type- 
I solution with N = n + 1/2 at L 2 = (3/4) {2n + l) 2 - m 2 /2 , where m = 
1, 2, 2n. For fixed values of k, n, m we find a single solution A 2 (k, (2n + 1) /m) 
obeying the following properties 

2 

Ttl 

A 2 (fc = 0, (2n + 1) /m) = 



6 (2n + l) 2 - 4m 2 
A 2 (A; = 1, (2n + 1) /m) = 1/4 (3.31) 

This is illustrated in figure 13.11 for two different values of the pair (n, m) by 
the solid curves, the dashed curve representing the limit (13.241) . 

The energy of the general type-Ill solution is given by the integral 



E 



in 



2m + 1 
V2A 



K(k) 



dy 



(ai - 1 + a 2 sn 2 (y, k)f + - (l + 16A 4 (k 2 - 1 - fc 4 )) 



which, by means of K(k), E(k) becomes 



(3.32) 
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Em 



2m +1 
V2A 




+ 




(3.33) 



For k = 1 one obtains 



E IU (l,m,n) 



4(2m+ 1) 
3v"2 



(1, m) . 



(3.34) 



This shows that solution III approaches solution II in the limit k — > 1, a fact 
which also appears in [88] . 

The energies of the solutions of some low-lying branches are presented in 
figure I3~2l as functions of L . For L > 5/2 all four types of solutions coexist 
and satisfy 



Ei (L, N = 0) < E I (L, N = 1) < E UI (L,m = l,n = 1) < E n (L, m = 1) < E (L) . 

(3.35) 



The four stars on the upper part of the figure show the position of the four 
bifurcation points L = 3.279, 3.775,4.093,4.272 of the n = 2, m = 3,2, 1,0 
solutions respectively from the N = 2 type-I solution. The two lower stars 
show the bifurcation values (L = 2.179, 2.5) of the n = 1, m = 2, 1 solutions 
respectively from the N — 1 type-I solution. 

3.5 Stability Analysis 
3.5.1 Type-I solutions 

In order to analyze the stability of type-I solutions we follow the steps done in 
[BE] which uses notions that can be found in |110j and |115j . Thus perturbing 
the fields <f) a , a = 1, 2 around the classical solution f l3.12l) -( !3~l~3l) . denoted here 



by 0: 



(p a (x) = 4> c a (x) + r) a (x) exp (-iut) 
we are led to the following equation for the normal modes: 



(3.36) 
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Figure 3.2: Goldstone model: The energies of some of the solutions are plotted 
as functions of the parameter L. The four stars on the upper branch indicate the 
four bifurcation points on the N = 2 type-I solution. The two lower stars indicate 
the two bifurcation points on the N = 1 type-I solution. The two numbers in 
parentheses refer to the number of negative and zero modes of the corresponding 
solution. 



A N,L 



>h 



to 



>h 



(3.37) 



A (N,L 



dx 1 



2 1 



where uj 2 is the eigenvalue and 




c = cos \Nx/L) , s = sin [Nx/L) , N + 1/2 = N 




(3.38) 



(3.39) 



We derived the above formula by substituting (13. 36ft into the field equations 
(13.81) . Terms 0(rj^) with s > 2 are considered small. 



The complete list of eigenvalues of the operator A yN, Lj for N > 1 
can be obtained by classifying its invariant subspaces. This can be done by 
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using the Fourier decomposition. Type-I solution has a twisted translational 
invariance and one can check that for an integer n > N the following finite- 
dimensional vector spaces 1 are preserved by A ( N, L 



( px px ( ~ \ 

V n = Span |a p cos — , (3 P sin — , p — n = ymod2Nj , \p\ < nj , (3.40) 



~ ( px px / ~ \ "] 

V n = Span |a p cos — , — j3 p sin — , p — n = \mod2NJ , \p\ < nj (3.41) 

under the following condition 

a p = p p if p - n + 2N > (3.42) 

where a p , f3 p are arbitrary constants. The operator A (jV, can then be 

diagonalized on each of the finite-dimensional vector spaces above, leading 
to a set of algebraic equations. 

To be more specific, define Ai = 2 (1 — A2), A2 = (n/L) . Also, consider 



cos ( ) x \ I cos ( ) x 



the vector 



(3.43) 

Acting with the operator A ^N, on one of the above vectors (say 
^jv+fe) we have the following steps: 



V^ +fc = 



r d 2 




2c 2 


2sc \ 


\ dx 2 




2sc 


2s 2 J 




r d 2 




1- A : 




\ dx 2 


2 



2 Ai 
2 



N+k 



(3.44) 



1 These vector spaces should have been the same as in (34) of [55] (with N — > iV) 
but they are not, due to a misprint in [55]. Here we correct this by writing these spaces 
explicitly in equations (|3.40j) and (|3.4ip . 
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Figure 3.3: Goldstone model: The values (|3.1ip are plotted as functions of L 
for j = 0,1,2 (solid curves), together with the values (1, 1, — 1), w|(l,0, — 1), 
w§(l,0, 1), (1,2,-1), w§(l,l,l), wg(l,3,-l) of (ET351) (dashed and dotted 
curves). The numbers indicate the multiplicity of the eigenvalues. 



where 1 is the 2 x 2 unit matrix and 



M 



cos O^jfj sin 
sin (f -cos 




which has the action 

MV„ +k = V„_ k (3.45) 

The same happens if we choose V^_ k instead. In that case, the equation 
above is MV^_ fc = V^ +k as expected. A matrix which has exactly the same 

Ai/2 



action as above can be found easily and this is 
Finally, we observe that 



Ai/2 
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which implies the condition that the determinant of the above matrix 
must be zero in order to acquire the non-zero eigenvalues. 

The eigenvalues of A ( N, L ) on V n are 



u 2 (N,k,±l) 



± 



1 

J? 
1 



N + 



+ k 2 



± 



\ 



L 2 




4k 2 [N + 



(3.46) 



for k = 0,1,2,... and similarly on V n for k = 1,2,3,.... Notice that if we 
replace N + 1/2 by N then the above result is the same with eq.(36) of [88] . 

For A" having a specific value and L slightly greater than N, the solution 
( 13. 121) - (13.1 31) possess 4N negative modes corresponding to k — 1, 2, 3, 2N 
in (13.461) (minus sign). When L increases, the number of positive eigenmodes 
increases as well 



L 2 = ^ (2N + l) 2 - m = 1, 2, 2N, (3.47) 

i.e. (cf( l3.27l) ) at those values of L where the solutions of type-Ill with n = N 
bifurcate from the solution of type-I. For 

L 2 > L 2 cr (iV) = 2 (2N + l) 2 - \ (3.48) 

all the modes are positive and (13. 121) - (13. 131) are classically stable solitons. 
These results are illustrated in figure [3T21 for N = 1 and N = 2. The num- 
bers in parentheses represent the number of negative and zero modes of the 
corresponding branch. 



3.5.2 Type-II and Type-Ill solutions 

The equations for the fluctuations about the solution (13.151) decouple to take 
the form of Lame equations: 

~ + 6k 2 sn 2 (y, fc) J % = n 2 lVl (3.49) 

~ + 2k 2 sn 2 (y, k)} V2 = n 2 2V2 (3.50) 
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where 



x = VT+¥y , n 2 a = (u 2 a + l) (k 2 + 1), a =1,2 (3.51) 

and uj 2 is the effective eigenvalue of the relevant operator. Equations (13.491) 
and (I3.50p admit two and one algebraic modes, respectively, with correspond- 
ing eigenvalues 

tt\ : 4 + k 2 , 1 + 4k 2 and Q 2 2 : 1 + k 2 (3.52) 

The corresponding values of uj 2 follow immediately from (13.511) : they have 
signature (+, +) and (0), respectively. 

It's a property of the Lame equation that the solutions determined al- 
gebraically correspond to the solutions of the lowest eigenvalues. The re- 
maining part of the spectrum therefore consists of positive eigenvalues. The 
spectrum of the equation ( 13. 49ft was studied perturbatively in |109j . while 
the relation between the Lame equation and the Manton-Samols sphalerons 
was first pointed out in [lll] - pT3j . 

Concerning the stability equation of type-Ill solution, it seems to be more 
difficult to deal with as the presence of the Jacobi Elliptic functions prevent 
us from having an analytic expression for Q(x) in (13.41) . Thus, we are not 
able to follow the steps done in the case of type-I solutions (i.e. find the M 
matrix) where trigonometric functions are easier to handle. Detailed analysis 
especially for the case of Lame equations can be found in |112] and |115] . 

3.6 Conclusions 

A detailed analytical study of the static solutions of the Goldstone model 
on a circle has been given in [SB] and we followed the same path here for 
our boundary conditions. Many results of [88] are connected with ours by a 
simple change on variables used. We write them explicitly above whenever 
is necessary. We also note our effort to impose mixed boundary conditions 
as well. Specifically, we enforced 0i to be antiperiodic and 2 periodic but 
there was no solution to satisfy this choice so it's needless to extend beyond 
this small remark here. 

Many details were presented on the stability analysis for the solutions 
found. Classically stable solitons were identified, together with the range of 
the parameter L for which they are stable. 
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This simpler model we present here and many others can be connected 
and can also give us the experience to deal with realistic (3 + l)-dimensional 
particle physics models in our search for possible metastable localized soli- 
tons. Ref. [88] on which this note was based, has connections with |109] and 
|110j as it concerns the branches of their solutions, also with [13] , while there 
are also interesting physical applications [114] as it is already mentioned in 
|88j . Soliton solutions in (3 + l)-dimensional models with antiperiodic bound- 
ary condition imposed on one spatial dimension, which is compactified on S , 
are analyzed in |115] - [TT7j where supersymmetry breaking is examined. 
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3.7 APPENDIX 



3.7.1 Part I: Jacobi elliptic functions 

The following serves as a small and helpful supplement on the mathematical 
functions we used above, as well as some of their basic properties, which were 
useful. 

Consider the following elliptic integral of first kind 

r do . 

u= / . (3.53) 
Jo VI - k sin 2 9 



Then we have the following definitions of the Jacobi elliptic functions 

sn u = simp (3.54) 
cnu = cos<p> (3.55) 



dnu = yl — A; sin ip (3.56) 

There are other nine Jacobi elliptic functions which are constructed from the 
above three but we don't use them here. The parameter k is called elliptic 
modulus and < k < 1, so the elliptic functions can be thought of as being 
given by two variables, the amplitude ip and the parameter k. 

From the above definition, it is easy to see that 

cn 2 + sn 2 = 1 (3.57) 

is valid. Also, 

dn 2 + k 2 sn 2 = 1 (3.58) 
The derivatives of these functions are 

— sn(z, k) = cn(z, k)dn(z, k) (3.59) 

— cn(z,k) = —sn(z,k)dn(z,k) (3.60) 

— dn(z, k) = —k 2 sn(z, k)cn(z, k) (3.61) 

where z, in fact, depends on u. 

Practically, sn(x, k) is the solution of the following differential equations 

+ (1 + k 2 )y - 2*Y = 0, ($L\ = (1 - y 2 )(l - k 2 y 2 ) (3.62) 
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and cn(x, k) is the solution of the following equations 

g + (1 - 2k")y + 2*Y = 0, =(l-y 2 )(l-e + kY) (3.63) 

while dn(x, k) is the solution of the following equations 

^f- 2^ + 2^ = 0, (%) =-(l- y 2 )(l-fc 2 -y 2 ) (3.64) 
dor \ ax I 
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4.1 Introduction 

The evolution of the Universe is believed to involve several symmetry break- 
ing phase transitions, out of which, topological defects, such as strings, are 
created [HI 11001 1106] . These transitions can be examined in the framework of 
condensed matter systems. Although there are differences from the cosmolog- 
ical case where, relativistic dynamics must be used and gravity is important, 
the formation of such defects in the laboratory [HH], can provide helpful hints 
for cosmology. Topologically stable knots and vortex-like structures in gen- 
eral, are of wide interest in condensed matter physics. For example, one can 
think of Bose- Einstein condensates (BEC) (i.e. see [8]), vortices in superfluid 
Helium-3 and Helium-4 [9] , or nematic liquid crystals [TTl [T2] . 

Also, in the framework of high energy physics, future experiments in 
LHC could answer whether metastable particle-like solitons exist in mini- 
mal supersymmetric Standard Model or two-Higgs Standard Model (2HSM) 
or not. In pJ)|,[T3], work on classically stable, metastable quasi-topological 
domain walls and strings in simple topologically trivial models, as well as 
in the 2HSM has been done. These solutions are local minima of the en- 
ergy functional and can quantum mechanically tunnel to the vacuum, not 
being protected by an absolutely conserved quantum number. One can also 
find other interesting subjects involving superconducting vortex rings such 
as, rotating superconducting rings [ID], electroweak strings [921 193] or work 
on such rings in SU(2) non-Abelian Yang-Mills-Higgs model [M]. Finally, 
twisted semilocal vortices examined in [95] can be connected to the models 
we present below while one can also search if stable rings can be made out 
of these strings. 

In this chapter, we consider a U(1)a model with a modified Ginzburg- 
Landau (GL) potential. The modification has to do with the addition of 
a cubic term. In Thermal Field Theory, such term comes from the 1-loop 
radiative corrections to the GL potential [102] . We search whether this model 
can admit stable strings or not. The features of such strings, if they exist, 
are the supercurrent which flows on the surface of the defect within a certain 
finite width, as well as, a magnetic flux in the interior of the defect (fig l4.ll) . 
This magnetic flux is a consequence of the existence of the supercurrent. 
One has to keep in mind that the magnetic field can penetrate in a certain 
depth inside the superconducting regions where the supercurrent flows. If the 
penetration depth is greater than the width of the superconducting surface, 
then the defect becomes unstable and can be destroyed. The GL potential we 
use here, is used in condensed matter physics as well (see [101] and references 
therein) . 
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Figure 4.1: U(1)a model: The relative position of the supercurrent as well as of 
the magnetic field on a xy-profile of the system on the right. The left picture is 
strongly reminiscent of an infinite solenoid. 

Stable defects of this U(1)a model, can also be used to form torus-like 
strings and study their stability. This can happen by taking a piece of such 
straight string and periodically connect its ends together. These string loops 
are examined in [SU] (chapter 5) but in the frame of a (7(1) x U(l) model, 
where the existence of the defect is ensured for topological reasons. That 
model, is a continuation of previous work [151 HB]- I n chapter 6, we will 
examine an extended version of the model in [90j (chapter 5). 



4.2 The U(1) A model 

This model consists of a complex scalar field i/j and a gauge field A^. The 
Lagrangian density describing our system is: 

C = -\f^ + \D^\ 2 -UM) (4.1) 

where the covariant derivative is = 8^ + ieA^, the strength of the field is 
Ffj, v = d^Ay — d u A^, while e is the U(1)a charge. The Lagrangian (14. ip . is 
invariant under the following U(l) gauge transformation 

^_ e »«^ A ^A,--d^) (4.2) 
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where 6(x) is a position-dependent phase. We choose the potential 

wi) = ^i^r^i^i 2 -fi^i+|j (4.3) 

where a, /3, 7 constants. We can set 7 = 1. Since y^y has dimensions of mass, 
we count the energy of the system in units of «/7- The vacuum is = 0. 
This vacuum leaves unbroken the gauge symmetry U(1)a- When 

w = w s ^Ei (4 . 4) 

(for /3 > 2), we have U(1)a 1 giving non-zero mass to A. Thus, one may 
generate an electric current flowing along regions where ^ 0. In fig j4.2l one 
can see the shape of the potential. The equation (14.41) gives the position of 
the minimum of interest for every (3 > 2. When /3 — 2, the secondary (non- 
trivial) minimum of the potential disappears, at = 1 position. When 
2 < (3 < = Pmt, h becomes zero only at = 0, while another minimum 
with non-zero \ip\ forms. When (3 = f3 cr it, the potential has another zero at 
I = \/2 which is also a local minimum. Finally, when j3 > {3 C rit it has two 
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more zeros at 



between which, it becomes negative (fig J4. 6ft . The mass spectrum is 



m A = 0, m v , = - (4.6) 



2 



4.3 The U(1)a model: Search for stable vortices 

We are interested in configurations having cylindrical symmetry, that is, 
infinite straight strings. The field ip can be non-vanishing on a cylindrical 
surface of specific radius. At infinity (p — * oo), we have the vacuum of the 
theory = 0. The ansatz for the fields is: 

j,(p, ip, z) = P{pY M \ A(p, <p, z) = ^A<p (4.7) 

P 

where M the winding number of the field ip and p, <p, z are the cylindrical unit 
vectors. We use cylindrical coordinates (t,p,tp,z), with space-time metric 
5V = diag(l, —1, — p 2 , —1). We work in the A = gauge. For the gauge 
field we suppose the above form based on the following thought: The A field is 
the one produced by the supercurrent flowing on the cylindrical surface. The 
current is in the (p direction thus, we expect the non-vanishing component 
to be and the amplitude P of ip to be independent of ip. As it concerns 
the scalar field ip, since it follows the geometry of the cylindrical defect, we 
expect that its amplitude is independent of z as well. 

With the above ansatz, the energy functional for minimization takes 
the form 



1 p2 



2 (dpA^Y + (d p P) 2 + — {eA v + M) + U(P) 



(4- 



E = 2vr / pdp 
Jo 

and the potential is 

U (P) = -P l I -P l — —P + - ) (4.9) 




The gauge field A has a magnetic field of the form 

Id A. 



B = (4.10) 
P dp 
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while, the field equations are 

V-- p • aP 

p p 

1 

p 

The usual rescaling arguments lead to the virial relation: 



1 P nP / \ 
d 2 p P + -d p P - -(eA v + M) 2 — [P 2 - PP + lj = (4.11) 

«9 2 ^ - -<9 P A, - 2eP 2 (eA v + M) = (4.12) 




27r / pdp\ 



2p 2 




Define 




2?r / pdp — - — - — + - 



For a solution of the model, we theoretically must have I1+I2 = 0. In fact, we 
define the index V = ffiljjlj ■ We want this index as small as possible. Other 
virial relations can be derived as follows. For example, one can consider the 
double rescaling of p — > Xp and either P — > /iP or — > or even both 
of the fields and then demand d\E\\ =1=fl = = d^E\\ =1=fl . 



4. 4 The £7(1)a model: Numerical results 

We use a standard minimization algorithm to minimize the energy functional 
( 14.81) . The algorithm is written in C. One can find details about the algo- 
rithm used, on page 425 of |104j but, briefly, the basic idea is this: Given 
an appropriate initial guess, there are several corrections to it, having as a 
criterion the minimization of the energy in every step. When the corrections 
at the value of the energy are smaller than w 1CT 8 the program stops and 
we get the final results. We are interested in final configurations having non- 
trivial energy. This signals the existence of a stable vortex with that energy. 
We check our results through virial relation (14.13j) . Finally, our results must 
also satisfy the field equations (14. lip . (14.121) . 

The initial guess we use for our computation is: 

P(p) = £iP M (i-tanh(0.2p 2 )) 

A v — - tanh(£ 2 p 2 ) 
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Figure 4.3: U(1)a model: In the first two graphs we have the initial guess (dashed 

lines ) as well as the final configuration of fields (solid lines — ) P(p), A v (p). 

We chose M = 1, e = 1, (5 = 2.17, a = 43.7. The energy E = 35.3 and virial is 
10 -4 . The bottom graph gathers all the fields. For the field P, the area between 
the horizontal lines P3 and P4 is energetically favorable (see also fig l4.6p . 



where £1,^2 are constants, the value of which, depends also on the location 
of the minimum (say = \ipo\) of the potential. For the final configurations 
we present in the figures l4~3l 14. 4[ we chose £1 = 1.75, £2 = 2. 

The initial guess also satisfies the appropriate asymptotics 

• near p = 0: P ~ p M , ~ p 2 

• at infinity: P — ► 0, exponentially 

while P must be non-zero somewhere between p = and p — > 00. One must 
make a careful choice of the initial guess. That is to say, the maximum value 
of P in the initial guess (dashed lines) must be inside the favorable area 
denoted by the horizontal lines in figs l4.3ir4~4l This area is dictated by the 
form of the potential and especially by its negative sectors (see fi.g J4.6l) . 
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Figure 4.4: U(1)a model: In the first two graphs we have the initial guess (dashed 

lines ) as well as the final configuration of fields (solid lines — ) P(p), A<p{p). 

We chose M = 1, e = 1, /3 = 2.13, a = 1104. The energy E = 115.2 and virial is 
2 • 10 -3 . The bottom graph gathers all the fields. For the field P, the area between 
the horizontal lines Pi and P2 is energetically favorable (see also fig !4.6p . 



Since, for values of (3 where 2 < (3 < (3 cr n is valid, we find no non- 
trivial solution, we searched and tried to find out what happens when we 
make the minimum of interest deeper. For that reason we searched in the 
region where (3 > flcrit- It is possible to find solutions to this model until f3 
reaches 2.13 (from above). Under this value, this is difficult if not impossible. 
Even at (3 — 2.13 (fjg J4.4[) we use great values of a in order to find the 
solution exhibited. The relation between a and (3 can be found in fi.g J4.7l 
The solutions of the model for two different values of (3 > f3 cr n, are shown 
in figs l4.3ll4~4l and a comparison between them in fi.g J4.5l in order to observe 
their different features. 
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Figure 4.5: U(1)a model: Here we have a plot in order to compare the changes 
of P(p) and Ba(p) as a goes from 2.17 (dashed and dotted - • - ), to 2.15 (dashed 
- - - ) and 2.13 (dotted • • •). 

4.4.1 Analysis for 2.13 < ft < 2.2 

We observe that for a specific ft, there is a small range for the parameter a of 
the potential, where the model exhibits the solution presented in the figures. 
Out of this small range and for greater a, we end up to a final configuration of 
negative energy. This happens because, the bigger the parameter a becomes, 
the stronger the potential is, thus the field P strongly prefers to acquire 
the value where the non-trivial minimum of the potential is (see fi.g J4.2p . in 
order to decrease further the energy. But, for ft > ft crit we have U < at 
the position of that minimum, which also enforces the total energy E to be 
negative in this case. We have to note that, virial relation of such a final 
configuration is not satisfied due to the great values the term (d p P) 2 acquires 
around p = and p — > oo. This is clear if one directly observes (14.131) . which 
can not be satisfied for U < 0. 

On the other hand, for smaller values of a, we get the trivial configura- 
tion (P = 0,^4^ = 0), as the benefit from the potential term is no longer 
satisfactory in order to have a non-trivial P. 

4.4.2 Reasons for instability when ft > 2.2 

Now, for high values of ft (i.e. ft > 2.2) we faced difficulties in finding a 
solution. We believe that this has the following explanation: as ft grows, the 
area of values where the potential U is negative, increases as well. 

If the parameter a is large, then the potential becomes a strong factor in 
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reducing the energy, as it is highly negative. Thus, for large a it is energet- 
ically favorable to decrease further the potential value. The latter happens 
by converting P function in such a way, so as to be, as much of it as possible, 
inside the energetically favorable area (it is denoted by the horizontal lines 
in the figs !4.3ll4~4~l) . Thus, we end up to a final configuration with potential 
U « and virial relations can not be satisfied (see for example eq J4.13l) as 
we have a sum of positive terms. 

On the other hand, for lower a the potential is no longer a strong factor 
for reducing the energy. In this case, it is energetically favorable to reduce 
the value of the (d p P) 2 and (dpA^) 2 terms as the system can gain more from 
these. The consequence is that P leaves the area of stability as its peak lowers 
in order to reduce the two terms above and there is only one possibility: to 
end up to zero energy, that is to say, the trivial configuration. 

The difference in the stable solutions we have found above, is that the 
values of (3 are such, that the potential is negative but not strongly negative 
while the changes on the terms (d p P) 2 and (d p A v ) 2 can deform the field P 
in such a way, so that it can still be inside the favorable area. Then, the 
potential term has the possibility to change in such a way, so it can satisfy 
virial as well. 



4.4.3 Reasons for instability when (3 < (3 crit 

In the following explanation we will use fig |4.5l & 14.61 In fig J4.5l one can 
observe that as we get closer to the critical value j3 = 3/\^2 = j3 cri t, the P 
field tends to acquire everywhere the value P = Pq (location of the non-trivial 
minimum of the potential). This leads to greater values of d p P. 

The above has a reasonable explanation which can be found in fig l4.6l 
As we get closer to /3 cr u, the space within the lines, where the potential can 
get negative values, becomes smaller. Under (3 crit the potential can be either 
positive or zero (the latter for P = only). 

Observe the energy functional to be minimized: 



E 



2/T / pdp 
'o 



P 2 



^(d p A v ) 2 + (d p P) 2 + — (eA v + M) 2 + U(P) 
zp p 



(4.14) 



The main target of minimization is to "fix" all the above terms in order 
to have the minimum possible value for the energy. Below, we analyze the 
possible cases. 



crif 
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Figure 4.6: U(1)a model: The P+ and P_ solution of eq. (|4.5|) . for which the 
potential becomes zero. Between these lines the potential gets negative values (see 
also figS2D. The plot is P vs. (3 for (3 > (3 crit . = -L. 



In that case, all terms except for the potential term, can be either pos- 
itive or zero. The potential term (as we saw in fjg J4.6[) can become 
negative for a range of values of P. But, as (3 decreases, this range be- 
comes narrower (as we see in fig |4.6l as well as in figs J4.3ll4~4l where this 
range is represented by the space between the two horizontal lines) and 
(d p P) 2 increases. This happens because the energy functional tends 
to decrease its value through the negative values of the potential term. 
We believe that this can not continue for (3 very close to f3 crit due to the 
fact that the range we described above, becomes so small, that P tends 
to get everywhere a constant value (the value Pq, which makes the po- 
tential negative). But P must be zero at p = and p — > oo, thus there 
will be a considerable increase in the (d p P) 2 term of the functional, 
and that makes the benefits of the negative value of the potential to go 
away, while the trivial solution P = becomes energetically favorable. 

• < Pa-it- In that case, the potential term can no longer become 
negative. It can be either positive or zero and because of the fact that 
the energy functional is a sum of five positive terms, it's reasonable to 
prefer the zero value which, at the same time, minimizes all the terms 
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Figure 4.7: U(1)a model: The plot exhibits the relation of the parameter a 
with respect to (3. We observe that as we approach the limit ficriu we need an 
increasingly deeper minimum of the potential which is expressed through the fastly 
increasing value of a. The dashed vertical line signals the position of p cr u = 

of the functional. 

From all the above, one can observe that the crucial difference between 
the above two cases of f3, is that in the j3 > f3 cr u case, the energy can have 
a minimum value (through the potential term) which corresponds to a non- 
trivial solution for P. The existence of negative values of the potential are 
the "way-through" that make this possible. 

4.5 Conclusions 

We studied a U(1)a model with a GL potential with a cubic term added 
to it. After the numerical analysis we did, we came to the conclusion that 
for (3 < Peru we find no non-trivial solution. For 2.13 < (3 < 2.20 we get 
non-trivial solutions which have the profiles we present in figs l4.31l4~4l Over 
(3 = 2.20 we have no non-trivial solutions. We analyze and explain our results 
in these cases. 

The form of the potential we have in this search, can be found in con- 
densed matter physics as well. On the other hand, one could try to make a 
loop out of the straight string studied above. We did this, but it was difficult 
to study mainly due to the fact that there are two instability modes, one 
having to do with the defect itself and another which has to do with the loop 
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that tends to shrink due to its tension. The former instability is excluded in 
[9U] for topological reasons. More details about the latter, can be found in 
the next chapter. 
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Chapter 5 

On axially symmetric solitons 
in Abelian-Higgs models 
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5.1 Introduction 

In a series of papers [THl EES] classically stable, metastable quasi-topological 
domain walls and strings in simple topologically trivial models, as well as in 
the two-Higgs Standard Model (2HSM) were studied. They are local minima 
of the energy functional and can quantum mechanically tunnel to the vacuum, 
not being protected by an absolutely conserved quantum number. In [JS] 
a search for spherically symmetric particle-like solitons in the 2HSM with 
a simplified Higgs potential was performed without success. Although the 
existence of spherically symmetric particle-like solitons in the 2HSM has not 
been ruled out, we shall here look, instead, for axially symmetric solutions 
in a similar system. 

Consider a model with superconducting strings [ISl |lU3j . Take a piece 
of such string, close it to form a donut-shaped loop and let current in it. A 
magnetic field due to the supercurrent will be passing through the hole of the 
donut (fig J5.2j) . The energy of the loop has, a term proportional to the length 
of the string and will tend to shrink the radius of the donut to zero and the 
ring to extinction. However, another force opposes this tendency. Namely, as 
the loop shrinks, the magnetic field lines are squeezed in the hole, since, due 
to the Meissner effect, they cannot leave the loop. They are trapped inside 
the hole of the donut, oppose further shrinking and might even stabilize the 
string. A brief plan of what we are searching, can be seen in fig J5.lL 

This, as well as other arguments [H]-[50] are inspiring but not conclusive. 
The magnetic field will not be trapped inside the loop if the penetration depth 
of the superconductor is larger than the thickness of the ring. Also, once the 
magnetic field gets strong it can destroy superconductivity and penetrate 
[105] . Further, there is a maximum current a superconductor can support 
(current quenching). This sets a limit on the magnetic field one can have 
through the loop, and this may not be enough to stabilize it. Thus, the 
above approach may work at best in a certain region of the parameter space, 
depending also on the defect characteristics [T5]. The purpose of this work 
is to apply the above straightforward idea to search for string loops in a 
U(l) x U(l) gauge model and to determine the parameter space, if any, for 
their existence and stability. 

Another interesting subject is to have a rotating ring. The rotation is 
another extra factor which could help the ring to stabilize. This work was 
done with success both analytically and numerically in (30] where vortons are 
exhibited. Another recent example of rotating superconducting electroweak 
strings can be found in [92], while for a review on electroweak strings, the 
reader should also check [93] . Finally, a work on static classical vortex rings 
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1 2 3 4 5 

Figure 5.1: U(1)a x C^(1)vk model: Step 1: Infinite straight string. Step 2: Cut 
a piece of it. Step 3: Periodically connect its ends together to form a loop. Step 
4A: Rings without current shrink due to their tension... Step 5 A: ...and finally 
collapse. Step 4B: Superconducting vortex rings can be stabilized?... Step 5B: ... 
or not? 

in SU(2) non-Abelian Yang-Mills-Higgs model can be found in [94"] . 

5.2 The U(1) A x U{l) w model 

The Lagrangian density of our model is: 

C = - \W% + \D^\ 2 + \D,<j>\ 2 - H) (5.1) 

where the covariant derivatives are D^tp = d^ip+ieA^ip, D^ip = d^+iqW^, 
the strength of the fields are F^ v = d il A v — d v A^ = d^Wv — d v W^ while 
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e and q stand as the relevant U(l) charges. The above Lagrangian (15. ip . is 
invariant under the following U(l) gauge transformations 



where b(x), p(x) are position dependent phases. We choose the potential U 

UM, M) = f (H 2 - .x 2 ) 2 + f (H 2 - vl) 2 + ||0|VI 2 - f 4 (5.2) 

The vacuum \(f)\ = v\ ^ 0, = 0, breaks U{l)w x ^ — >■ ^ (1)a, giving 
non-zero mass to W. The photon field stays massless. There, U(v\,0) = 0. 
The vacuum manifold A4 in this theory is a circle S 1 and the first homotopy 
group of M. is iz\{M.) = iii(S 1 ) = Z ^ 1 which signals the existence of 
strings. In regions where |0| = 0, the field is arranged to be non-vanishing 
and U(l)w x U(1)a — > U{l)w Thus, U(1)a — * 1 and one may generate an 
electric current flowing along regions with vanishing |0|. Hence, this theory 
has superconducting strings |103j . 

In fact, with \<p\ = F and = P, the extrema of the potential are four: 

• F = = P with U(0, 0) = \g x v\. 

. F = 0, P = v 2 with 17(0, v 2 ) = \g x v\ - \g 2 v\. 

• F = v 1 ,P = with U{v u 0) = 0. 

• F = a\, P = a>2 with U(ai, a 2 ) = — ^g^f 



with 



a 3 



\ 1/2 / \ 1/2 

, g2(g\v{ - gsvl) \ I gi(&vl - g 3 vf) 

d\ — g 5 °2 — 2 

^1^2 - gi J \ g\g2 - gi 

gz{g\gz{gz - 9\9^\ + g^gg - g\g2)v\ - ^gig^jgj - g\g2)v\v 

Kgm - gl) 2 



The vacuum that is of interest to us is: |0| = v\, \ip\ — and leaves unbroken 
the electromagnetic U(1)a- 
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We can choose <p = (/i(x) + v^e 101 ^' and ip = / 2 (x)e l62 ^ with < 
0|/i(x)|0 >= =< 0|/ 2 ( x )|0 >. Plugging these in the Lagrangian, we 
have the masses of the fields: 

m A = 0, m w = qvi, = g x v\, = -(fl^i - g 2 vf) (5.3) 

When |0| = v\, the potential becomes: 

U(v u m = \i;\ 2 (^\^ + k^j (5.4) 

where k 2 = ^(fl^i — 92V 2 )- Searching for the minimum we find: 

dU 



0^ \i>\ \2k 2 + g 2 \ij\ 2 =0 (5.5) 



When k 2 > 0, we have only one minimum at = 0. Otherwise (fc 2 < 0) we 
have a minimum at |^| 2 = — which means that the photon acquires mass. 
Thus, we choose the condition k 2 > which reads as the right equation that 
follows: 

v I) 2 

9i > #2-f, #3 > 92^ (5.6) 

The left equation is another condition which comes from the fact that when 
|0| = and | ip\ = f 2 the potential has value 

= (5.7) 

and in order to ensure that this is greater from the value of the potential at 
the vacuum, which is 0, we arrive at the condition above. This makes the 
value of the vacuum of interest, lower than the other extremum values as 
well since 

£7(0,0) > U(0,v 2 ) > = C/(0,0) (5.8) 



5.3 The U(1)a x U(l)w model: Search for superconducting vor- 
tex rings 

Configurations with torus-like shape, representing a piece of a U(l)w — * 1 
Nielsen- Olesen string |18j . closed to form a loop, are of interest in this search. 
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Thus, we will require <p to vanish on a circle of radius a (the torus radius) 
<f)(p = a, z = 0) = 0. At infinity (p — > oo, z — ► oo), we have the vacuum of 
the theory. This translates to |0| — > Vi, \ip\ — > 0. The ansatz for the fields 



is: 



<Kp,^,z) = F(p,z)e me ^ 
il>{p,<p,z) = P(p,z)e m * 

MP,<P, Z ) = — 

P 

W(p,<p,*) = W^(p^) z 

where M, N are the winding numbers of the relevant fields, p, tp, z are the 
cylindrical unit vectors and we define 

Q(p,z) = arctan ( J (5.9) 



We use cylindrical coordinates (t,p,<p,z), with space-time metric g^ u = 
diag(l, — 1, — p 2 , — 1). We work in the A = = W° gauge. Especially 
for the gauge fields, we suppose the above form based on the following rea- 
sonable thoughts: The W-field is the one related to the formation of the 
string thus, it exists in the constant-<p plane. This means that in general its 
non- vanishing components will be W p and W z . The A-field is the one pro- 
duced by the supercurrent flowing inside the toroidal object. The current is 
in the p direction thus, we in general expect the non- vanishing component to 
be A,p. Finally, as it concerns the scalar fields 0, ip, they follow the geometry 
of the toroidal defect which has axial symmetry thus, we expect that their 
amplitude is independent of p. 

A more general choice for would be 

0(p, p, z) = F(p, z yMe( P ,z)+i x ( P ,z) ( 5 _ 1Q) 

where x(p, z), an arbitrary function. But gauge invariance allows us to 
change (p — > (pe lb ^ l \ where 6(x) an arbitrary space-dependent phase. We 
can choose fo(x) = —x(p, z ) thus, gauge fixing removes the arbitrary function 
X- 

We proceed by extracting the energy- momentum tensor and by consider- 
ing the above ansatz we will derive the energy functional. Couple the system 
to an external gravitational field g pv and define 

= £=4^- ( 5 - n ) 

v -g dg^u 
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r 5 



z 



Figure 5.2: U(1)a x U{V)w model: A x — y profile of the superconducting ring 
(left) as well as a y — z profile (right) where one can view how the mechanism we 
propose, against shrinking, works. 



The model coupled to the background gravitational field is 
£ = 



- \g^g^F, u F Xp + g^(D^)*D^ - \g^g vp W^W Xp 



and by using 



we get 



T ap = ^F^g^ + 2{D a ^)*D^-W^W^g^ + 2{D a < 1 



f-9 



(5.12) 



We are interested in the energy of the system which is the integral of the 00 
component of the energy-momentum tensor 



E = J d n x^Tg T 00 



(5.13) 
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where yfg is the determinant of the spatial part of the metric. The energy 
density for the static ansatz above is 

T oo = lg*gfl FijFkl _ ^il), c yn jC + IgVWijWu - gV&WDrf + U 

(5.14) 

The terms of the energy density are 

KVFF 1 ( {d p A,f {d z A,f 



-(d p W z - d z W p ) 



P 2 



-^(D^yDjij = (d p Py + (d z Py + —(eA v , + N) 2 
-g^(Di4>)*Dj(f> = (d p F) 2 + (d z F) 2 + ((qW p + Md p Q) 2 + (qW z + Md z Q) 2 )F 2 
Thus, with the above ansatz, the energy functional takes the form: 

POO POO -j . 

E = 2n Vl / pdp / dz —((d p A v ) 2 + (d z A v [ 
+ (d p P) 2 + {d z P) 2 + P ^{eA v + N) 2 + 



(d p F) 2 + (d z F) 2 + -(d p W z - d z W p f + 
{qW p + Md p Q) 2 + {qW z + Md z O) 2 ^F 2 + U(F, P) 
and the potential U (fig J5.3p can be written as follows: 



U(F,P) = |(F 2 - l) 2 + |(P 2 - u 2 ) 2 + |F 2 P 2 - |« 4 



(5.15) 



(5.16) 



where u = v^jvy. This is the energy functional we use in our numerical work. 
The conditions ( 15.61) to be satisfied by the parameters become: 

9i > 92U 4 , g 3 > g2U 2 (5.17) 
The gauge fields A, W have magnetic fields of the following form: 

l (dAp , dAp 

V x A = B A = - I z p 



p\ dp 



dz 



V x W = B 



w 



dW z _ dW p 
dp dz 



s6 
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The field equations follow: 



d l F + d i F+ 9 f> F 9ii 



P 



1 )F-—P 2 F 



qW p -M 



zcos 2 Q 



qW z + M 



cos 2 
(p-a) 



2i 



F = 



d n P 



-\ 2 P 92fp 2 _ u 2^p_93_ F 2p = Q 



eA^ + N, 2 



dlW p - d z d p W z - 2qF 2 \qWp-M 



2eP 2 [ eA v + N 
zcos 2 Q 



1 



(P ~ «) 2 
cos 2 Q 








We can also write down the currents associated with field, namely j£ 
and jf and the total current out of these as well as the supercurrent 
associated with the ip field. These are 



T 



i 2p P 2 

= V (it) 2 + tit ) 2 , X* = ——( & A V + N) (5.18) 



where 



-2qF 2 (qW p + Md p B), jf = -2qF 2 (qW z + Md z Q) (5.19) 
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Figure 5.4: U(1)a X U(1)w model: For different values of y/gi, we plot the lower 
bound over which the condition (I5.20p is satisfied. The plot is P max vs. e and we 
plot the function P ma x 



ILL 

e 



As explained in the Introduction, Meissner effect is of crucial importance 
for the stability of the torus-like object. The magnetic fields produced by the 
supercurrent penetrate into the toroidal defect in a distance dictated by 
the penetration depth. In this theory, the C/(1)a symmetry breaks inside the 
string and the photon acquires mass m\ = e 2 < P > 2 where < P > is the 
expectation value of the charge condensate in the vicinity of the string core. 
There is the superconducting sector of the defect. The penetration depth is 
A = — = /L . But < P >= Pmax < u thus, one can have an estimate for 

mA e<P> mux — j 

A, with a lower bound for its value. This is A > — where the equality holds 
when Pmax = u. On the other hand, one can also compute an upper bound 
for the thickness of the defect. We know that r& = — = —=. If our concern 
is to search for stable rings, a reasonable step is to demand the penetration 
depth to be smaller than the string thickness which means 



\<r (j> 



eP 



< 



2 7 }2 . 



(5.20) 



This is the condition needed in this case. From the above condition, we 
get the diagram shown in fig |5. 41 where one can see the area where it's more 
possible to find stable solutions if they exist. The numerical results we present 
later, are what we found while searching inside this region. 
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-F l3 F l3 + |A0| 2 + -rWijWij + |AVf + U{\<t>\, = £ (5.21] 



Below, we analyze a way to derive some virial relations in order to check 
our results. The energy density for a static solution, which is of interest here, 
reads: 

\f^ + \dm 2 + \ 

where = diAj — djAi and Wij = diWj — djW{. For the solution of the 
system, we can write: 

5E 5E 5E 5E 5E 5E , 

» = W = wrw = W = JZ =0 (5 - 22) 

Define 

rp X rp xrp xrp tp 

which together with 

SE dT 00 n dT 00 . , 

can be written in a shorter form /j = o^Cr^ where 

9s: 9s de 

G - = £5i >-WW) di *~W^) d4 *~W^) diAk 



0(c^) 9(9,-^*) 9(9,-Wz] 
This means that any static solution of the field equations satisfies 

djGij = (5.25) 



But, we have 



de 
d{dj(j>*) 

de 
d(djA k ) 

de 

de 
de 



Drf 

djAk — dkAj = Fjk = eijkBAi 
Drf 

djWi - diWj = Wji = e ifl B w . 
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which leads to 



Ga = eSij-iDrfYdtf-DjWF-WjAWi- 

- (Djipydrt - Djipdiip* - F jk d t A k = 

= sSij - (DjWDtf - + iqW4{D^)* - 

- iqWifDrf - WjiWa - WjAWi - 

- (DjipyDiip - Dj^D^y + leA^D^)* - 

- ieAiip*Djip - F jk F ik - F jk d k Ai = 

= eb i3 - (bi4>ybj(j) - {brfybrf - w u w jt - ^(w^wi) 

+ Wt (dtWji - iq((f)*b j( j) - - 

- {D^yDjiP - (D^yD^ - F jk F ik - d k (F jk Ai) + 
+ A { (d k F jk - ieWDtf - (D^yi;) 



where in the second equality we just add and subtract the necessary terms in 
order to have the terms only with covariant derivatives. The latter happens 
because all the other terms are gone as we see below: 

d J d k (d j A k -d k A,)A i = => (dp k A k -d 2 k d j A j )A i = 
d j d l {d j W l -d l W j )Wi = U (d]d l W l -dfd j W J )W i = 



and from the field equations, the following terms are also zero 

5E 5T 00 

= = -iqiPDrf - ftDrf)*) + diW jt = 



6Wj d{d l W j ) 
SE 5T 00 

thus 

Gi^eSy-WuWji-ibityDrf-iDrfybi^^ 

(5.26) 

but 

Fi k Fj k = e nik B A„e m jkB A m = ^A^ij ~ Ba^a^ 
WuWji = e p uB Wv e q jiB Wq = B^5y — B W .B W . 
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which means that 

Gij = ( - S _ 3l + \D(f>\ 2 + |Zty| 2 + Uj + B Wi B w . + B Ai B Aj - 
- (Diip)*Djip - (Dj-ip)* Di<ip - (D i( j))*D j( f> - {D j( j))*b4 

From (15.251) one obtains 

J d 3 xG ik = = J d'xdjixkG^) = J dS jXk G i3 (5.27) 

The last equality follows from Gauss theorem, while the first is due to the 
fact that 

dx k 

dj^XkGij) = -Q^r G ij + x k (djGij) = 5 kj Gij = G ik (5.28) 
where we also used (15.251) . Take the trace of (15. 2 7ft . We have 

J d 3 xTrG tk = 2vr J pdpdz ( - S _ Si + \D(f)\ 2 + \D^\ 2 + 3EM (5.29) 

while the left hand side of ( 15. 27ft in the (p direction of integration of the 
toroidal object gives 

2tt J dS v G w = (5.30) 
and (15.271) finally ends up to the following: 

2tt J pdpdz [B 2 W + Bl)-i?J=0 (5.31) 



Another way to derive virial relations is through Derrick's scaling ar- 
gument. The virial relation for the field of our model, must have the 
constraint 4> K (p = a,z) = <fi{p = a, kz) = 0. Consider the rescalings p — > p, 
z — > kz, F K — > F, P K — > P, A Vk — > W p , Zk — > kW P;Z . Then, we find the 
minimum through the relation W- = when k — 1. The virial relation for 

° OK 
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our model, which is a way to check our numerical results is: 



poo POO 

Ii = 27TV\ / pdp / dz 

Jo J-oo 



\{d P W z - 8 Z W P ) 2 + 7^(d z A v ) 2 + (d z P) 2 



+ {8 z Ff + 2F 2 qW p {qW p + Md p Q) + {qW z + Md e G)' 



POO POO 

-2nvi / pdp / dz 

JO J -co 



1 

Yp 



2 (d p A,) 2 + (d p F) 2 + (d p P) 2 



+ ^{eA v + N) 2 + {d z W p ){d p W z -d z W p ) + 



+ F 2 (qW p + Md p G) 2 + {qW z + Md z G) 2 + 



i±( F 2 _ 1)2 + 92 { p2 _ n 2 )2 + 93 F 2p2 _ 92 4 



and we must have h+F = 0. We define the index V = ^K 1 ! , l 1 ^) and we want 

I-11I+I-/2I 

its value to be as small as possible. We can derive many other virial relations 
by assuming generally for a field <f>, the "double" rescaling 4>(x) —>■ n<j)(px) 
and then demanding ^| re=1=jU = = ^\ K= \ =p . For example, consider the 
following rescalings p — > p, z — > pz, F K — > F, P K — > «P, A, 



.4 



- W p , 2 . We have 



+ 



/*oo /*oo 

27rui / pdp / dz 

Jo J -oo 



2p- 



.{d z A^) 2 + {d z P) 2 + {d z F) 2 + 



+ 2F 2 Md z G{qW z + Md z G) + -(0„W, - d z W p 



2f 



POO POO 

—2ttvi / pdp / 

JO J~oo 

^{eA v + iV) 2 + d p W z (d p W z - d z W p ) + 



(<VM 2 + {d p Pf + (<9 P F) 2 + 



+ F 2 [(qW p + Md p Q) 2 + (qW z + Md z Q) 2 ) + 



+ 9_ 1{F 2_ 1) 2 92 {p 2_ u 2 } 2 + i± F 2 p 2_92 * 

4 v ' 4 K 1 2 4 
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together with 







poo 


POO 


h 


= Attvi 


\ pdp 


/ dz 








— oo 






POO 


POO 


h 


= 2txv\ 


/ pdp 


I dz 






h » 


— oo 



(d p P) 2 + (d z P) 2 + — (eA v + Nf 



^2 D 2 



where, as above, we must have I3 + I4 = = J5 + I§. 



5.4 The U(1)a x U(l)w model: Numerical results 

A standard minimization algorithm is used to minimize the energy functional 
of f)5.15p . The algorithm is written in C. One can find details about the 
algorithm on page 425 of |lU4j but, briefly, the basic idea is this: Given 
an appropriate initial guess, there are several corrections to it, having as a 
criterion the minimization of the energy in every step. When the corrections 
on the value of the energy are smaller than ~ 10~ 8 the program stops and we 
get the final results. A 90x20 grid for every of the five functions is used, that 
is, 90 points on p-axis and 20 on z. Here, we begin with fixed torus radius 
a. Then, the configuration with minimum energy for this a is found. Other 
values of a are chosen as well and the same process goes on until we plot 
the energy vs. the torus radius E(a). It would be very interesting to find a 
non-trivial minimum of the energy (in a m « n 7^ 0), which would correspond to 
stable toroidal defects with radius a min . 

The initial guess (figure 15. 5j) we use for our computation is: 

tanh((p-a) 2 + ^ 2 ) M / 2 
tanh(p 7V )(l - tanh((p - a) 2 + z 2 ) 

~7 tanh {((p-ly + z^) 

Mz cos 2 ef (p-a) 2 + z 2 \ 
q(p-a) 2 y (p - a) 2 + z 2 + (a 2 /4) ) 

_Mcos 2 ©/ {p-af + z 2 \ 
~ OTP -<>■) \(P- ^ 2 + z 2 + (a 2 /4) J 

where £ a constant. Figures [5.61 15.71 indicate what exactly we plot in fig J5.5l 
and where. This initial guess also satisfies the appropriate asymptotics 



P(p,z) = 

Ap(p,z) = 

W p (p,z) = 

W z (p,z) = 
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Figure 5.5: U(1)a X U(1)w model: A typical plot of the initial guess we use for 
the five fields, for the lowest winding state M=l,Af = lonz = plane. 




Figure 5.6: U(1)a X U(1)w model: This is a typical picture of the toroidal soliton 
we are after. 
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Figure 5.7: U(1)a x model: The profile of the torus has two identical 

cycles. The initial configuration of fields (plotted in fig J5.5[) is the same in both of 
them thus, we choose to plot only the right side. 



near p = 0: 



F^O, P~p l \ A v ~p 2 f(z) 



near (p = a, z = 0): 



F~p M/2 , W p = = W z 



(5.32) 



(5.33) 



at infinity: 



l-0( e -Vp) j p^C(e 
M 

9p6 oo + e 

9 



M 
9 



a) 2 + z 2 . 



where p = (p 

For fixed torus radius (i.e. here a ~ 2.2) we present a typical graph of 
the final configuration of the lowest energy (see fi.g J5.8l) . We also present the 
plot of the energy of the system vs. the radius of the toroidal object which 
reveals the instability of the system as there is no non-trivial minimum (see 

figEHD. 

From the results, we can point out a few things. Firstly, the greater 
the value of e we use, the stronger the supercurrent becomes. Secondly, the 
greater the value of e we use, the lower the radius a where the supercurrent 
quenches (fi.g J5.9p . These are expected as the increase of e makes the con- 
dition of equation (15.201) stronger, something which means that the mass of 
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1 2 3 4 5 6 



P 

Figure 5.8: U (I) a xU(l)w model: Typical plot of the final configuration of fields. 
Parameters are M = 1, N = 1, e = 5, q = mw = 2, g\ = 14, = y/gi = 3.74, 
g 2 = 12, g 3 = 14, u = 1, = y/(g 3 - ff2« 2 )/2 = 1, v x = 7.5 • 1(T 3 . Energy 
E = 0.72 and virial is 1.5 • 1CT 3 . 



the photon increases and the penetration depth decreases at the same time. 
It is also reasonable that a stronger current can "defend" the defect, against 
the magnetic field, a little longer. 

Another observation is that as g 2 increases and becomes close to g±, we 
see that the maximum current increases and the quenching takes place again 
at lower a (fig J5.10p . This is expected as one can see from the potential in 
equation (15.1 6p of the energy functional, because the stronger the coupling g 2 
is, the more important the relevant term g 2 {P 2 — m 2 ) 2 /4 becomes. The latter 
has as a consequence, the increase of P max which counteracts the effects from 
the increasing g 2 coupling. 

The parameter space where we searched, starts from g\ = 4 (m^ = 2). 
In order to search the model, we reached values around g\ = 30 (m^ ~ 5.48) 
over which, e has to be very large in order to respect the condition f)5.20p . 
There is also the fact that great values of e in general is an unwanted feature 
since we use a semiclassical approach. As it concerns the other couplings we 
have g-3 = g\ and g 2 = g\ — k with 0.5 < k < 8, (u — 1). 
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Figure 5.9: U(1)a x U(l)w model: The top graph is the energy vs. the torus 
radius E(a). The middle graph is the supercurrent 1^ vs. a. In that graph one 
can clearly see current quenching. The bottom graph is the quantity e 2 P^ iax 
vs. a (or rriA vs. o) where one can see the area in which the condition of 
eq J5.20l holds (lines over the g\ -limit line). The resistance from the magnetic 
field can be seen as a sharp increase on the supercurrent. Dotted lines are for 
e = 6, dashed for e = 8 while dashed and dotted for e = 10. All plots are for 
(M,N,u, 9l ,g 2 ,g 3 )=(l, 1, 1, 14, 12.5, 14). 



5.5 Explanation concerning the instability of the vortex ring 

From the condition of equation (I5.20p . it is understood that we are enforced 
to lower gi as much as possible and/or increase e. But these steps are not 
as easy as they might seem. There are some limitations. The lower bound 
on the value of gi has a reasonable explanation. For low values of g±, the 
coupling g 2 is also low (because g\ > g2U A ). Now, when g 2 is small enough, 
the changes on the term g 2 [P 2 — u 2 ) 2 /4 are unimportant for the energy, in 
comparison to the term (d p P) 2 . In that case, the lowering of the last term 
minimizes the energy, something which means that P — > 0. Indeed, this is 
numerically observed. There is also a lower limit on e which is reasonable 
because the lowering of e results to an increasing penetration depth. 
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Figure 5.10: U(1)a x U(l)w model: The top graph is the energy vs. the torus 
radius E(a). The middle graph is the supercurrent vs. a. In that graph one 
can clearly see current quenching. The bottom graph is the quantity e 2 P^ iax vs. 
a (or wia vs. a) where one can see the area in which the condition of eq l5,20l 
holds (lines over the gi-limit line). The resistance from the magnetic field can 
be seen as a sharp increase on the supercurrent. Dotted lines are for gi = 12, 
dashed for §2 = 12.5 while dashed and dotted for cj2 = 13. All plots are for 
(M,N,u,e, gi ,g 3 )=(l, 1,1,6, 14,14). 



We searched for stable objects for values over these limits described 
above. The numerical results are exhibited in figs l5.91lo.lll and as we see, 
these objects are unstable. We argue that the explanation for the instability 
is current quenching and that, only for high values of e. For lower values 
of e, that is, of the order of 1, we have, according to equation (15.201) . that the 
penetration depth is much bigger than the string thickness thus, stability is 
out of the question. The latter is also numerically observed. 

Now, we base our aspect about quenching on qualitative as well as quan- 
titative arguments. We observe that as the torus shrinks, the supercurrent 
suddenly drops to zero which signals the destruction of the defect. Just be- 
fore the sharp drop, we notice that the supercurrent rises with increasing 
rate. This must be due to the resistance the torus meets from the magnetic 
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Figure 5.11: U(1)a x U(1)w model: The plot of the energy of the system 
vs. the radius of the torus for four different sets of parameters. Dotted is 
for (ffi, <72, 53)= (14,12.5,14), dashed is for (18,15,18), dashed and dotted is for 
(25, 20, 25) while dashed with two dots is for (30, 24, 30). For all data sets we have 
(e, M, N, u)=(10, 1, 1, 1). As one can observe, there exists no minimum. 

lines as it shrinks. One can observe that as the supercurrent increases and the 
condition of equation (I5.20p is satisfied at the same time (i.e. see dashed and 
dotted line of fig J5.9p . suddendly the current is lost. This can be explained 
only through current quenching. The above phenomenon is not observed 
when equation (15. 201) is not satisfied (i.e. see dotted line of fig J5. 10ft . There, 
as the magnetic lines penetrate the ring, they meet almost no resistance since 
A is much greater than r^. 

Another observation which supports our quenching argument is that, 
rough estimations on the maximum current a string can sustain, lead us to 
the following formula (see page 129 of [105] or Appendix I at the end of the 
chapter) which makes a small overestimation in order to have an upper limit: 

It ax < vW (5.34) 

where a = J J dpdz P 2 . In figs l5. 9115.101 the maximum value of the 
supercurrent is close to the limit of the estimation of equation (I5.34p . The 
table below gathers the estimated lf sL (according to equation (15.340 ) and the 
computed maximum supercurrent (Zf om ) for the parameters of these figures. 
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Finally, one can make an estimation of the value of the supercurrent 
which would stabilize the ring, namely 1^ tab ■ This can be done as follows. 
As explained in the introduction, there is the tension of the string which 
shrinks the loop and the magnetic field which opposes this tendency. When 
the ring is stabilized we have E tension = E magnetic . Here, E tension ~ 1 and 
Emagnetic = 2irvi f pdp J dz(B\/2). Without any calculation, one can point 
out that, since the total energy in the quenching radius is below E tension = 1, 
then the E magnetic which is a fraction of it, would be even smaller. Recall 
that Ba oc Z^, which means that we need a Tf tab which is well above the 
maximum current we can have inside the defect. Calculations of the mag- 
netic energy are in agreement with the above observation and place its value 
around E magne u c ~ 5 • 10~ 3 << E ten sion ~ 1- This translates to the following 
conclusion: lf tab > 10 -lf sL . 

Thus, we find out that the current needed for stabilization, is at least ten 
times bigger than the value of the maximum current our string can sustain. 
We also observe that our numerical maximum current values are close to the 
theoretical estimations about quenching. This means, that we will have cur- 
rent quenching as an "obstacle" towards stabilization, since the supercurrent 
will not be enough in order to create the magnetic field needed. 
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radius of the torus a 



Figure 5.12: U(l)w model: Profiles for the scalar field F(p,z) on the left and 
for the magnetic field By/ on the right side, for three random values of a. On 
the bottom graph there is the energy vs. radius of the torus o. Parameters are 
(M,m^,mvy,fi)=(l, 4.47, 2,0.01). Black solid line is for z = 0. 

5.6 U(l)w model: Vortex ring without supercurrent 



Here, we present a simpler model of a string of toroidal geometry but without 
the ip scalar field responsible for the existence of supercurrent. In this model, 
we only have the scalar field 4>(p, z) together with its gauge field W(p, z). The 
presence of the latter gives us the freedom to choose Q(p, z) = arctan(z/(p — 
a)), as we did before. The initial guess we use for the fields is the following: 



W p {p, z) 
W z (p,z) 



tanh((p-a) 2 + z 2 ) M / 2 
Mz cos 2 ef (p-a) 2 + z 2 



q(p-a) 2 \(p-a) 2 + z 2 + (a 2 /4) 
ios 2 9 / (p - a) 2 + z 2 



q(p -a) \(p- a) 2 + z 2 + (a 2 /4) 



As there is no current flowing inside the string there will be no Meissner 
effect which means that there is no reason to prevent the vortex ring from 
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collapsing. Thus, we expect to find no non-trivial minimum on the energy of 
the system, something which signals the instability of that object (fig J5.12l) . 

5.T Conclusions 

Future experiments in LHC could answer whether metastable particle-like 
solitons exist in MSSM or 2HSM or not. In [15] there is a search for spher- 
ically symmetric solitons in the frame of the 2HSM with a simplified po- 
tential. Here we search for axially symmetric solitons which, if stable, will 
have a mass of the order of TeV [T5]. We considered the U{1)a x U(1)w 
model, where the existence of the vortex is ensured, for topological reasons. 
There, we searched for stable toroidal strings. We present and analyze our 
observations. This chapter tries to answer to expectations having to do with 
observations of stable axially symmetric solitons which would be possible to 
detect in later experiments of LHC. For relatively small values of e (~ 1), 
which are of interest in that case, the system seems to have no stable vortex 
rings. In fact, this instability is present even in other parameter areas where 
we searched (i.e. e > 6 see figs l5. 9115. 101) . We explain why we believe that 
the main reason of instability is current quenching. 
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5.8 APPENDIX 



5.8.1 Part I: Derivation of the formula T^ xax < y/aeu 

To derive the above formula, we follow the steps of |105] (pages 129-130). 
The current density is given by 



f = ie(ipD^ - ipD^ip) (5.35) 

where ip = P(p)e 1 *^ and the covariant derivative = <9 M + ieA^. Sub- 
stituting above, the absolute value of the current density is j v = 2eP 2 D v & 
and the total supercurrent is the integral of the cross section of the current 
density 

/*oo poo 

1= dp dz f (5.36) 

JO J -co 

from which we get I = 2aeD ip Q, where we denote J J dpdzP 2 by a. 

The term g2{P 2 — u 2 ) 2 /4 of the potential, enforces the expectation value 
of P to reach u. On the other hand, in the energy functional there is the 
term P 2 (D V §) 2 which tries to "counteract" the above potential term. The 
latter means that (_D^$) 2 acts as a "negative mass" squared. Because of that 
term, the expectation value of P, say Pq, is not u but 

P 2 « M 2 - — (Z^$) 2 (5.37) 

92 

The above decrease has an effect on the supercurrent as well. It decreases 
its value thus, 

(Z^$) 2 



1 = 2ae(D lp $) 



1 - 2a- 



v? 



(5.38) 



In the above equation, only D v § has dimensions (of mass) and we rename 
it as X. We want to find the maximum value of the current: 

dX u , 



and by substituting, we find out that the maximum value of the supercurrent 
is 

Imax < V&eu (5.40) 
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5.8.2 Part II: Toroidal coordinate system 

The toroidal coordinates (77, £, <p) are related to the cartesian (x, y, z) through 
the following definition: 



a a a 

x = — smh 7) cos ip, y = — smh 77 sm <p, z = — sin £ 



D 



D 



D 



(5.41) 



with a a scale factor of the torus, D = cosh 77 — cos£ and < 77 < 00, 
< £ < 2tt, < if < 2tt. 

In order to find the field equations, we need to relate the derivatives of 
cartesian coordinates with those of toroidal coordinates, thus the following 
formulas can be very useful: 



d cos <yj(cosh?7 cos £ — 1) d sin £ sinh 77 cos tp d Dsiwp d 

dx a dr) a <9£ a sinh 77 dip 

d siny?(cosh77cos£ — 1) d sin £ sinh 77 sin tp d Dcosip d 

dy a dr) a <9£ a sinh 77 dip 

d sin £ sinh 77 d cosh 77 cos £ — 1 <9 

dz a dr] a 



d* + d 2 + d 2 = -— 



(cos £ — cosh 77) 



<9 2 d 2 
+ 



sinh 77 <9y? 2 <9£ 2 



+ 



+ sm ^-^z + (cos 4 cosh 77 - lj-r-; — tt- 
a£ smn 77 ar/ 



c) 2 



<9t7 2 



The relations between unitary vectors of the toroidal and the cartesian 
coordinate system follow: 



h 



D 



cosy?(l — cosh 77 cos £) i + sin<^(l — cosh 77 cos £) j — sinh 77 sin £ k 
1 



D L 



sinh 77 sin £ (cos </? i + sin tp j) + (cos £ cosh 77 — 1) k 

ip = — sin ip i + cos 99 j 



5.8.3 Part III: Cylindrical coordinate system 

The cylindrical coordinates (p, tp, z) are related to the cartesian (x, y, z) 
through the following definition: 



x = pcostp, y = psmtp, z = z 



(5.42) 
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In order to find the field equations, we need to relate the derivatives of 
cartesian coordinates with those of cylindrical coordinates thus, the following 
formulas can be very useful: 



d d snap d 



dx dp p dcp 

d . d cos 9? d 

dy ^ dp p dip 

sin lip d 2 sin 2 <p d sin 2 cp d 2 

p dpdip p dp p 2 dip 2 

sin 2ip d 2 cos 2 ip d cos 2 <p> d 2 

p dpdip p dp p 2 dip 2 

cos 2<p d 2 sin 2 tp d sin 2 <p d 2 

p dpdip Ip dp Ip 2 dip 2 

v 2 - -—( —\ 1 92 

p dp \ dp J p 2 dip 2 dz 2 

The relations between unitary vectors of the cylindrical and the cartesian 
coordinate system follow: 

p = cos ip i + sin ip j 
ip = — sin tp i + cos ip j 
z = k 



d 2 


d 2 


sin2y9 d 


dx 2 




p 2 dip 


d 2 


d 2 


sin2(/) d 


dy 2 




p 2 dip 


d 2 


sin2ip d 2 


cos2<yj d 


dxdy 


2 dp 2 


p 2 dip 



5.8.4 Part IV: Field profiles in different z levels 

Here, we present the dependence of the fields of the U(l) x U(l) model 
examined above, on z variable. For that reason, we plot the fields of fig J5.8l 
for the same parameters, but in different z levels in order to observe their 
behavior. 
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Figure 5.13: U(1)a x U(1)w model: Typical plot of the final configuration of 
fields. Solid line: z = 0, dotted line: z = 0.6, dashed line: z = 1.2, dashed and 
dotted line: z = 1.8. Parameters are M = 1, N = 1, e = 5, g = = 2, 
5i = 14, = y^I = 3.74, g 2 = 12, g 3 = 14, u = 1, m^, = ^(#3 ~ 92U 2 )/2 = 1, 
Wl = 7.5 • 10~ 3 . Energy E = 0.72 and virial is 1.5 • 10~ 3 . 
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Chapter 6 

The extended £7(1) x U(l) model 
for vortex rings 
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6.1 Introduction 

Stable strings of the previous U(l) model (examined in chapter 4) could be 
useful in order to create vortex rings and study their stability. It could also be 
relatively helpful numerically, as we would have four fields for minimization 
in the energy functional instead of five we had in [90] (chapter 5). This is not 
as easy as it might seem, since there are two instability modes. The vortex 
itself is not necessarily stable while forming a torus, and the latter has the 
tendency to shrink due to its tension. The first instability can be avoided in a 
U(l) x U(l) model as the one presented in [90J (chapter 5), where a numerical 
search for bosonic superconducting static vortex rings in a U{1)a x U{1)w 
model is examined. There, the existence of straight strings is ensured for 
topological reasons. The superconductivity of the loop though, does not 
seem to prevent shrinking. The conclusion there, is that current quenching 
takes place before stabilization. 

Here we deal with an extended version of the previously discussed model 
of chapter 5. We add appropriate higher derivative terms which might help 
the ring to stabilize. We present and analyze our results. 

6.2 The U(1) A x U(l) w model 

The model discussed in [HO] (chapter 5), is being described by the Lagrangian 
density 

A) = -\f% - \wl„ + \D^\ 2 + \D,<f>\ 2 - M) (6.1) 

where the covariant derivatives are D^ip = d^ip+ieA^ip, D^(p = d^+iqW^, 
the strength of the fields are = d^A u — d v A^ = d^W u — d u W^ while 
e and q stand as the relevant U(l) charges. The potential U is 

U(\4>\, H) = |(H 2 - vlf + f (H 2 - vlf + f |0| 2 H 2 - |, 2 4 (6.2) 

The vacuum of this theory is \<p\ = v\ ^ 0, \ip\ = and breaks U(l)w x 
U (1)a — > U (I) a, giving non-zero mass to W. The photon field stays massless. 
There, U(v\,0) = 0. The vacuum manifold A4 in this theory is a circle S 1 
and the first homotopy group of M. is ixi(M) = 7ri(5' 1 ) = Z which signals 
the existence of strings. In regions where \<j>\ = 0, the field is arranged 
to be non- vanishing and U(l)w x U(1)a —> U(l)w- Thus, U{1)a 1 and 
electric current flows along regions with vanishing \<p\. Hence, this theory has 
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superconducting strings |103| . The vacuum of the theory leaves unbroken 
the electromagnetic U(1) A . For g 3 v 2 > g^v\ this vacuum is stable, while 
g\v\ > g2v\ ensures that it is the global minimum of the potential. The mass 
spectrum is 

m A = 0, m w = qvi, = g x v\, m\ = -(g^i ~ 92vf) (6.3) 
6.3 The extended U(1) A x U(l) w model 

We intend to modify the above model by adding higher derivative terms of the 
fields <p and ip and find out whether such changes can stabilize the ring. By 
following Derrick's scaling argument [98], one can argue that terms such as 
l-D^I 4 or l-D/i^l 4 or \D ^%Ij\ 2 could be helpful. Also, in an investigation 
of a similar model [97], the conclusions lead to the same path, in order to 
search for possibilities of stabilizing such solitons against radial shrinking. 
Thus, we have the following Lagrangian density: 

C = Co + c^\D^\ A + c^D^l* + c^\D^\ 2 \D^\ 2 (6.4) 

where c^, c^, constants. This Lagrangian (16. 4p with the extra terms, 
exhibits the symmetries of the original Lagrangian £ (16.11) . 

Configurations with torus-like shape, representing a piece of a U(l)w — > 
1 Nielsen-Olesen string, closed to form a loop, are of interest in this search. 
Thus, we will require to vanish on a circle of radius a (the torus radius) 
0(p = a, z = 0) = 0. At infinity (p — > oo, z — > oo), we have the vacuum of 
the theory. This translates to —* V\, — > 0. The ansatz for the fields 
is: 

<t>(ft,<p,z) 
ip(p,ip,z) 

A(p,ip,z) 

W{p,(p,z) 

where M, N are the winding numbers of the relevant fields, p, <p, z are the 
cylindrical unit vectors and the phase 

9(p, z) = arctan ( ) (6.5) 

\p — a) 



= F(p,z)e iM@ip ' z) 
= P(p,z)e iN ^ 

P 

= WJp,z) p + W z (p,z) z 
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We use cylindrical coordinates (t,p,tp,z), with space-time metric that has 
the form g in ,=diag(l, —1, — p 2 , —1). We work in the A = = W° gauge. We 
follow the ansatz of [901. 



With the above ansatz, the energy functional to be minimized takes 
the form: 

poo poo -f 1 

E = 2<kv x pdp dz —Ud p A ip ) 2 + (d z A^ 2 )+-(d p W z -d z W p ) 2 + 

JO J-OO V / A 

+ (d p P) 2 + (d z P) 2 + (d p F) 2 + (d z F) 2 + ^-(eA v + N) 2 + 

+ [{qW p + Md p Q) 2 + ( ? W 2 + M<9 Z 9) 2 ) F 2 + 

+ c^|(9 p F) 2 + (^F) 2 + ((giy p + Ma p 6) 2 + (giy, + M^e) 2 )F 2 | + 

+ c^(d p P) 2 + (d z P) 2 + y(eA 9 + iV) 2 | + 

+ c^j ^(d p F) 2 + (d z F) 2 + (few, + M<9 P 6) 2 + feW* + M^6) 2 )F 2 j 



(<9 p P) 2 + («9 Z P) 2 + — (eA v + N) 2 )\+ U(F, P) 



(6.6) 



and the potential U follows: 



U(F,P) = ^(F 2 -lY + ^(P 2 - u 2 Y + |P 2 P 2 - ^u 4 (6.7) 



.2 , 92 



where u = v-ijv\. This is the energy functional we use for our numerical 
analysis. The conditions to be satisfied by the parameters become: 



9i > 92U 4 , g 3 > 92U 2 



(6.8) 



The magnetic fields are 



VxA = B, 



1 dA M „ dA u 



p\ dp 



V x W = B 



w 



dW z _ dW p 
dp dz 
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while the currents associated with <fi field, namely j% and jf and the total 
current X^ out of these as well as the supercurrent X^ associated with the ip 
field are 



2eP 2 



(eA v + TV) 



(6.9) 



where 



^ = -2gF 2 (giy p + M9 p e), = ~2qF 2 {qW z + M<9 Z 9) (6.10) 



Finally, in order to check our numerical results, we can derive virial 
relations through Derrick's scaling argument. Below we present the virial 
relations we use in our search. Consider the rescalings p — > p, z — > kz, 
X K ► X, P K > P, -> Ay, W p ^ -> By demanding §f = when 

k = 1 and if we define 



I\ = 2tiv\ J pdp I dz 



\{d p W z - d z W p f + ^ 2 {d z A v f + (d z P) 2 



+ {d z F) 2 + 2F 2 \qW p {qW p + Md p Q) + {qW z - 

+ c^|4X 2 (^X) 2 [{qW p + M«9 p 6) 2 + + Md z Q) 2 J + 

+ 3(d z F) 4 + 2{d p F) 2 {d z F) 2 + 

+ ( 4X 4 ( + M<9 P 6) 2 + (qW z + M<9,6) 2 



• (W P (gI¥ p + Md p e) + (gW, + Md z Q) 2 yj + 
+ 4X 2 ((6> P X) 2 + («9,X) 2 ) (qW p (qW p + M<9 p 9) + (qW z + Md z Qf 



cA 3(d z P) 4 + 2{d p P) 2 {d z P) 2 + ^{d z P) 2 {eA v + N) 2 \ ^ 
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+ c^l {d z P) 2 {d p Ff + (d p P) 2 (d z F) 2 + 2>{d z P) 2 {d z F) 2 + 



p2 ( 2P 2 F 2 

+ ^(e^ + A0 2 (^F) 2 +f— — 



(eA v + N) 2 + 2F 2 ((d p P) 2 + (d z P)' 



qW p (qW p + Md p Q) + (qW z + Md z Q) + 



+ (d z P) 2 F 2 (qW p + Md p Q) 2 + (qW z + Md z Q)' 



POO f 

L pdp i- 


oo 

dz 




oo 



1 

2p 



2 {d p A,) 2 + {d p Ff + (d p P) 2 + 



P 2 
I 1 



{eA v + N) 2 + {d z w p ){d p w z - d z w p ) + 



+ F 2 {qW p + Md p Qf + {qW z + Md z ®f 



+ cJ {d p Ff + F 4 |^ p + Md p Q) 2 + (gW 2 + Md z Qf 
+ 2F 2 ((<9 p F) 2 + (^F) 2 ) [{qW p + M«9 P 6) 2 + + Md z Qf 



+ c^<( (d p P) 4 + *—(eA„ + iV) 4 + ^-(d p P) 2 (eA^ + N) 2 } + 



pi 



+ c^<( («9 p P) 2 (<9 p F) 2 + — (eA., + iV) 2 (9 p F) 2 + 



p2p2 



— (eA v + N) 2 + F 2 (d p P) 2 



{qW p + M<9 p 9) 2 + {qW z + M<9,9)' 



4 v 7 4 v 7 2 4 



we must have 7i + Ii — 0. We define the index V = an d we want its 

value to be as small as possible. We can have many other virial relations by 
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assuming generally for a field 0, the "double" rescaling <p{x) — > n<j)(fix) and 
then demand ^\ K=1=fl = = ^\ K=1=fM . For example, we check our results 
through the following relations as well. Consider the following rescalings 
p ^ p, z -> pz, F K -> F, P K ^ ftP, A Vk A^, W PiZ(t W Pi2 . We define 



J 3 = 27ri>i ^ pdp / dz 



— (d z A ip ) 2 + (d z P) 2 + (d z Ff + 



+ cJ3(d z FY + 2(d p F) 2 (d z F) 2 + 



+ 4F 4 ( ( 9 W P + Md p Q) 2 + {qW z + Md z Q) 2 j (Md z Q(qW z + M<9 Z 6)) + 



+ AF 2 {d z F) 2 {qW p + Md p Qf + (gW* + M<9 Z 6) 2 + 



+ 4F 2 ((d p F) 2 + («9,F) 2 ) (Md z ®{qW z + M<9 Z 9) 



3(d z P) 4 + 2{d p P) 2 {d z P) 2 + ^!(9 z P) 2 (e^ + Nf U 



+ C<H> { {d p P) 2 {d z Ff + {d z P) 2 {d p Ff + 3(^P) 2 (^P) 2 + 



P 



+ — (e^ + iV) 2 (<W + 



p2 

P 2 



(d z P) 2 + —(eA ip + N) 2 j(A 



+ (d z P) 2 F 2 (qW p + Md p Q) 2 + (qW z + Md z Q)' 
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POO f 

L 


oo 

dz 




oo 



yi^f + {d p Pf + (d p F) 2 + 



+ ^(eA v + N) 2 + d p W z (d p W z -d z W p ) + 
+ F 2 I (qW p + M<9 p 6) 2 + (qW z + M<9 Z 6) 2 ) + 



+ cA (dpFy + F (qW p + MdpQY + (lW z + Md z e)' 



+ 2F 2 ((o> p F) 2 + («9,F) 2 ) ^(g^ p + M«9 P 6) 2 + (qW z + Md z Gf 

{ P 4 IP 2 1 

+ cJ {8 p Pf + — (eA v + iV) 4 + — (^P) 2 (e^ + iV) 2 U 

+ c J (o> p P) 2 (<9 p F) 2 + ^(eA„ + iV) 2 (9 p P) 2 + 



p2p2 



{eA v + Nf + (d p P) 2 F 2 



{qW p + Md p e) 2 + {qW z + M<9,9)' 
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together with 



poo poo p2 

h = 4^/ pdp dz (d p P) 2 + (d z P) 2 + —(eA v + N) 2 + 

JO J-oo P 

f2P 4 
2{d p Pf + 2{d z Pf + —{eA, + Nf + A{d p P)\d z P) 2 



+ 



4p2 
„2 



(9 p P) 2 + (9 z P) 2 )(e^ + iV) 2 l + 



+ CH) \ (d p P) 2 (d p F) 2 + (d p P) 2 (d z F) 2 + {d z P)\d p Ff + {d z P) 2 {d z Ff + 



P 2 



P 2 ( (<9 P P) 2 + (d z P) 2 + —{eA v + JV)' 



{qW p + M<9 p 0) 2 + + Md z Gf 
+ ^(e^ + iV) 2 ((9 p P) 2 + (9 z P) 2 )| 
/ 6 = 2 ™if pdp J dz g 2 (p 2 -u 2 ^P 2 + g 3 F 2 P 



where, as above, we must have I3 + I4 = = 1 5 + 1$. 

In the energy functional (16.61) . the terms that come from the |P;0| 4 extra 
term, are multiplied with c^. In fact, these terms are proportional to d p F 
and d z F. Thus, if one chooses (c^, c^, c^) = (1, 0, 0), then there are more F- 
derivative terms in the functional. Energy minimization lowers these terms, 
something which one expects to lead to a thicker string. This is a wanted 
feature in order to stabilize the ring. This would have another consequence. 
The extra "P-terms", enforce the P field to stay away from its vacuum 
expectation value within a larger area. This means that P ~ inside a 
bigger area. There, the potential becomes 



(6.11! 



and its value increases because of the term gi(F 2 — l) 2 /4, as P — > 0. Mini- 
mization tends to make P — > u, which tries to compensate for that increase. 
The latter means that P will increase and this is another wanted feature 
since the supercurrent oc P 2 . 
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On the other hand, in (15.151) . the terms that come from the |Pj^| 4 extra 
term, are multiplied with c^. In fact, these terms are proportional to d p P 
and d z P. Thus, if one chooses (c^, c^,c^) = (0, 1,0), then the P-derivative 
terms become more. The minimization of them is expected to decrease the 
charge condensate P. This decrease is an unwanted feature. 

Finally, the terms that come from the \Di(j)\ 2 \Di%l)\ 2 term, are multiplied 
with c^. The consequences of the addition of this term can be seen if we 
observe that the extra terms are of the form 



{dpy 



{d p Ff + (d z F) 2 + 



(6.12) 



where the dots represent the rest of the terms which are positive as well. 
We expect a stronger decrease of the charge condensate P. This is because 
dP < 1 which means that (<9P) 2 > (<9P) 4 . Thus, the need for minimizing 
the derivative terms of P becomes stronger than in the case of (c^,, c^,, c^) = 
(0, 1, 0). Apart from this, it is also the fact that, since dF ~ 1, the weight of 
the P-derivative terms is now greater than unity and this is another factor 
which would tend to make P — > or, at least, smaller than in the case 
(c^c^c^) = (0, 1,0). 

Theoretically, the fact that the terms |Pj^| 4 and l-D^I 2 ^^! 2 tend to 
shrink the charge condensate P, can also be seen from the virial relation 
-^5 + Iq — above. The integral Iq is negative (since P < u) and the addition 
of extra terms leaves it unchanged. On the other hand, 1$ is a sum of positive 
terms and the above two extra derivative terms rise the value of 1$. Thus, 
the only way to satisfy that virial relation is either to increase g% and/or to 
decrease P. If J5 is large enough, then P will be enforced to become zero in 
order to satisfy the I 5 + J 6 = relation. 



6.4 The extended U{1)a x U{1)w model: Numerical results 

We use the same minimization algorithm, as in [5D], to minimize the energy 
functional presented in (16. 6p . A 90 x 20 grid for every of the five functions 
is used, that is, 90 points on p-axis and 20 on z. We begin with fixed torus 
radius a. Then, the configuration with minimum energy for this a is found. 
Other values of a are chosen as well and the same process goes on until we 
plot the energy vs. the torus radius E(a). It would be very interesting to find 
a non-trivial minimum of the energy (in a m i n 7^ 0), which would correspond 
to stable toroidal defects with radius a m j„. One crucial check of our results 
is done through virial relations. 
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Figure 6.1: Extended U(1)a X ?^(l)w model: A typical plot of the initial guess 
we use, for the lowest winding state M = 1, N = 1 on z = plane. 

The initial guess (fig JbMj) we use for our computation is: 

F(p,z) = tanh((p - a) 2 + z 2 ) M/2 
P{p, z) = tanh(p 7V )(l - tanh((p - a) 2 + z 2 ) 

a , ^ N u ( P 2 



Mz cos 2 Q (p-af + z 2 



q(p - a) 2 \ (p - a) 2 + z 2 + (a 2 /4) 



M cos 2 e 



(p - a) 2 + 2 2 



g(p-a) \ v (p-a) 2 + ^ 2 + (a 2 /4) 
This initial guess also satisfies the appropriate asymptotics 
• near p = 0: 



F^O, P~p JV , A^p 2 f{z) 



near (p = a, z — 0): 



F~p M/2 , W„ = = W* 



(6.13) 
(6.14) 
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F(p,0) 



1 2 3 4 5 6 




1 2 3 4 5 6 



Figure 6.2: Extended U(1)a x U(1)w model: A typical graph which ex- 
hibits the effects of the higher derivative terms on the scalar fields. Solid 
line is for (c^, c^, c<^)=(0, 0, 0), dotted line for (c^, c^, c^)=(0, 1, 0), dashed 
for (c^, C0, C0^)=(1, 0, 0) while dashed and dotted for (c^, c^, c^)=(l, 1, 0). In 
the case (c^, c^, c^)=(0, 0, 1), P is trivial. Parameters in this figure are 
(51, 92, 93, e, q, u, v u M, N)= (14, 12, 14, 6, 2, 1, 7.5 • 1(T 3 , 1, 1). 

• at infinity: 

F ~ 1 - C(e-^), P ~ (^(e^v 7 ^) 

Wp^-^eu + oce-^), ^--^eu + ^e-^) (6.15) 

where p = (p — a) 2 + z 2 . 

Based on [90] (chapter 5), we search on parameter areas where e acquires 
relatively large values, but they are interesting concerning the possible stabil- 
ity of the loop. The reason was analyzed in that paper (chapter) and stems 
from the need to have string thickness greater than the penetration depth as 
well as strong supercurrent. 

The results confirm our expectations stated previously. For example, 
in fig l6.2l one can observe that when the extra higher derivative term is 
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Figure 6.3: Extended U(1)a X f/(l)w model: The effect of the higher 
derivative terms on the energy and the supercurrent when these are plot- 
ted vs. the radius of the torus. Solid line is for (c^, c^, c t /^ ) )=(0, 0, 0), dot- 
ted line for (c^, c^, c^)=(0, 1, 0), dashed for (c^, c^, c^)=(l, 0, 0) while dashed 
and dotted for (c^, c^, c^)=(l, 1, 0). For (c^, c^, c<^)=(0, 0, 1), the supercur- 
rent as well as the charge condensate, are trivial. Parameters in this figure are 
(gi,g 2 ,93,e,q,u,v 1 ,M,N)=(U, 12.5, 14, 10, 2, 1, 7.5 • 1(T 3 , 1, 1). 



\Diip\ A , then F exhibits no change while the charge condensate P decreases 
(dotted line in fig J6.2l) . In this case, the consequence is the reduction of the 
supercurrent (dotted line in fig J6.3p when compared to the case of the original 
model without extra terms ((c0,c^,c^) = (0,0,0), see solid line in fig |6.3l) . 
On the other hand, when the extra term is \Di<f)\ 4 , then F widens and this 
leads to the broadening of P as well. The latter increases (compare dashed 
and solid lines of fig J6.2p and the supercurrent increases too (dashed line of 
fig |6.3l) . Finally, when the extra term is |-Dj0| 2 |-Dj^| 2 , then P = for the 
values of #2 we use. In general, the charge condensate can be non-trivial for 
higher gi. This happens due to the term g2{P 2 — u 2 ) 2 /4 of the potential. 
When #2 grows, P tends to reach u. This is also numerically observed. 

One can combine two extra terms to see what happens. In example, 
we add both |-D;0| 4 and l-D^I 4 ) (case (c^, c^, c^) = (l, 1, 0) in fig J6. 21) . This 



119 



The extended U(l) x U(l) model for vortex rings 



results to the addition of the "favorable" F-derivative terms, but also P- 
derivative terms would be present. This translates to the growth of P but 
not as much as in the case (c^, c^, c^) = (l, 0, 0). We also observed that the 
combination of either |-Dj0| 4 or l-Dj^l 4 or even both, with the |A0| 2 |A^| 2 
term, leads to shrinking of P because of the strong effect of the last term. 
After the theoretical and numerical analysis, we conclude that the most "in- 
teresting" extra term is |.Dj</>| 4 . 

It is clear that current quenching is present here as well (fig J6.3p . The ex- 
tra term |Dj0| 4 , can increase the supercurrent and can make the penetration 
of the magnetic field more difficult, as the string increases its diameter, but 
this increase is not enough in order for the ring to stabilize. For the shake 
of research, we also tried higher values of in order to make the favorable 
term more significant. We also tried higher values of e, but the ring could 
not stabilize in a non-trivial radius. 



6.5 Discussion 



The most crucial terms of the energy functional which could provide for the 
stability of the ring in a non-zero radius a, are: 




B = / pdp dz —(eA v + N) 2 

Jo J-oo P 

These two terms have an explicit total 1/p behavior which helps them to 
increase as the torus radius decreases. The problem is that they are not 
increasing at a satisfactory rate in order to overcome all the rest terms of 
the energy which are decreasing with p. This is also numerically observed. 
Under some radius a, the ideal would be to have a strongly increasing charge 
condensate P. Then, as the radius decreases and at the same time P in- 
creases, the above two terms would start to increase with sufficient rate in 
order to lift the energy of the system. The magnetic term would increase 
because as P 2 /p 2 increases, \A V \ — > \A Vma J — > N/e. The latter would make 
the derivatives of (and B& as well) to increase as p decreases. 

This is what we tried to do here, especially with the help of the |-D;0| 4 
extra term. The charge condensate became more robust but that was not 
enough. This supports the conclusion of [90] which states that very high 
values of supercurrent are needed for stabilization. It seems that such highly 
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increasing currents can not be produced, despite the help of extra terms. 
In fact, numerical details reveal that the rate of increase of the terms A, B 
above, is A(A + B) ~ 1CT 3 , while the rate of decrease of all the rest terms is 
A(E — (A+B)) ~ 3- 1CT 2 . This means that the rate of increase of A, B should 
be ~ 30 times bigger. That case would require P > u which is something 
that does not seem to satisfy I5 + Iq = virial relation. But even if that was 
possible, current quenching would be another "obstacle" . 

6.6 Conclusions 

We are based on the model of [90] (the previous chapter) and analyze an 
extended version of it, by adding higher derivative terms in order to check 
whether they can stabilize the superconducting ring or not. Although the 
|.Dj0| 4 term is helpful on that direction, it turns out to be insufficient and 
current quenching prevails. Finally, we discuss what one would need for a 
stable ring. This discussion in combination with the results of [90], seems to 
exclude the possibility of existence of such vortex rings in this model. 
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T.l Conclusions 

Bosonic superconducting cosmic strings and solitons were the main subject 
of research in this Thesis. Firstly, we examined a Goldstone model for an- 
tiperiodic solitons on S 1 . An analytical study as well as stability analysis 
are performed. Classically stable solitons were identified. Such model can 
be connected and can give us the experience to deal with realistic particle 
physics models in our search for possible metastable localized solitons. 

Then, a detailed numerical search for bosonic superconducting straight 
strings in a U{1) model with Ginzburg-Landau potential with a cubic term 
added to it, is done. Such strings exist in a small, numerically determined 
region. We fully analyze and explain the reasons of stability there, as well 
as the reasons of instability in the rest of the parameter space. Such models 
can be found in Condensed Matter Physics as well. On the other hand, 
this model can provide the basis for searching torus-shaped solitons in the 
framework of High- Energy Physics. 

The stable strings of the above U(l) model though, are not necessarily 
stable while forming a torus. Thus, we worked in a U(l) x U(l) model which 
admits straight bosonic superconducting strings due to topological reasons. 
We cut a piece out of these strings in order to create a torus and examine 
its stability. We performed a detailed numerical study which showed that 
such solitons seem to be unstable. An explanation concerning this instability 
is given with the help of theoretical predictions about current quenching as 
well. There are strong indications that current quenching is the main reason 
of instability. 

After the above conclusion, we searched further the model by extending 
it. We add higher derivative terms which could be in favor of stability of this 
vortex ring. We analyze and explain our expectations. Despite the fact that 
the results are improved, we again have non-stable objects. We discuss the 
reason for this and explain what it would be necessary for stability. Out of 
this discussion, the conclusion which states that these objects are unstable, is 
strongly supported. Current quenching seems to be the obstacle, while there 
are indications that the vortex rings of this theory, need some specific extreme 
conditions in order to be stable. Conditions which can not be satisfied by 
the U(l) xU(l) model we examined, even if current quenching phenomenon 
was absent. 

LHC experiments, which are about to start in early 2008, can answer 
whether metastable particle-like solitons exist in MSSM or 2HSM. In the 
framework of 2HSM, the existence of spherically symmetric solitons was ex- 
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amined but the result was negative. The next step was to search for the 
possibility of the existence of axially symmetric solitons with mass of the 
order of TeV. The confirmation or rejection of the existence of such solitons, 
would be a very interesting, important as well as urgent task, since LHC will 
start working in a few months in a very promising energy range. The work 
appeared in this Thesis, gives an answer and completes the above mentioned 
search on solitons with axial symmetry. 

Apart from this, search for such solitons can be done in other models with 
higher symmetries or extended models [H2]-[SI]- Experimental condensed 
matter physics is another branch where these gauge theories can apply. Some 
examples can be found in [3]-[l2] and [1001 HOlj . On the other hand, such 
theories can give detailed and useful hints for cosmological observation of 
strings in general and of loops in particular [2Zj-[SI], [HI]-[En|, [S3], |118| - [T2"0"] . 
Finally, another interest comes from the recent developments in Superstring 
theory. In the framework of large extra dimensions, long superstrings may 
be stable and appear at the same energy scale as GUT scale cosmic strings. 
More can be found in [ED E2], [125]- [153] . 
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NUMERICAL METHODS 

In this part, we present details of the numerical methods used in models 
of Part I. Then, we present full details of the algorithm we finally chose 
for our research in Part II. We tested the algorithms on models with known 
solutions such as, the Nielsen- Olesen vortex and the straight superconducting 
string. Newton- Raphson and energy minimization were comparatively easier 
to handle, while the latter was faster and quite helpful as it concerns the 
physics of the model and for that reason we chose it for our research. 
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8 . 1 Relaxation 



• How it works 

In relaxation methods, ODE's (Ordinary Differential Equations) are replaced 
by approximate FDE's (Finite-Difference Equations) on a grid or mesh of 
points that fills the domain of interest. The general rule is: 



where h = (x max — x m i n )/N and N stands for the number of points used to 
divide the region we are interested in. In other words it is the step on x-axis. 
The information the algorithm needs as input in order to start computing the 
final configuration of fields are, the initial and final conditions as well as an 
initial guess for the solution. It is important to have an, as good as possible, 
initial guess. This can be done by using the results of an asymptotical analysis 
on the behavior of the differential equations around x m i n and x max . An 
iterative process is used and the initial guess improves step by step, while 
every next step uses the improved configuration produced by the previous 
one. Finally, the result is said to relax to the true solution. 

• Concluding remarks for Relaxation 

The satisfactory output for Nielsen-Olesen was created after many difficulties 
since the algorithm was very sensitive on very small changes of the boundary 
conditions. This had effect on the model of superconducting bosonic string 
where the method did not led us to a satisfactory result. In general this 
algorithm was not easy to use and it's not efficient |104j especially in problems 
which have solutions with oscillatory parts. We continue by working with 
Runge-Kutta. 

8.2 Runge-Kutta 

• How it works 

This method advances a solution from x n to x n+ i = x n + h through the for- 
mula Un+i = Un+h f(x n , y n ). User supplies the code with the initial boundary 
conditions only. This means that the values of y Xmin and its derivative at that 




(8.1) 
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point are necessary. The "core" of the code of a fourth-order Runge-Kutta 
method follows: 

ki = h f(x n ,y n ) (8.2) 

/ h ki\ 
k 2 = h f [x n + -, y n + — J (8.3) 

h = hf(x n + (8.4) 
k 4 = h f(x n + h,y n + k 3 ) (8.5) 

y n+ i = y n + k i+ k i + k i + k i + o^) (8.6) 

o 3 3 o 

which means that every step h involves four evaluations of the right hand side 
as one can see. Apart from that, we also used the fifth-order Runge-Kutta 
method in order to achieve greater accuracy although our problems involve 
smooth functions, thus one could support that it might not be necessary. 
The " core" of the fifth-order Runge-Kutta method is presented below: 



h = hf(x n ,y n ) (8.7) 

k 2 = h f(x n + a 2 h,y n + b 2 ik 1 ) (8.8) 

k 3 = h f(x n + a 3 h, y n + b 3 ih + b 32 k 2 ) (8.9) 

k 4 = h f(x n + a A h, y n + &41&1 + ^2^2 + b A3 k 3 ) (8.10) 

h = h f(x n + a 5 h, y n + b 51 k ± + b 52 k 2 + b 53 k 3 + 6 54 /c 4 ) (8.11) 



k 6 = h f(x n + a 6 h,y n + b 61 ki + b 62 k 2 + b m k 3 + b m k 4 + b 65 k 5 \8.12) 
y n +i = y n + C]k\ + c 2 k 2 + c 3 k 3 + c A k A + c 5 k 5 + c 6 k 6 + 0(h e ) (8.13) 

where a^, bij, Ci are the Cash-Karp parameters and can be found on page 
717 of [TO!] . 

• Concluding remarks for Runge-Kutta 

Runge-Kutta was easier to use in comparison to Relaxation but the results 
were poor in the "test" -models we tried to use it. Thus, we proceed on 
Newton-Raphson and energy minimization algorithms. 



128 



8.3. Newton-Raphson 



8.3 Newton-Raphson 

• How it works 

The algorithm that follows, was created in 2004, in collaboration with Dr. 
E . D . M . Kavoussanaki . 

Suppose a differential equation written in the form of matrices: A-x = B, 
where A is a N x N matrix, x is a N x 1 matrix of the unknown function 
and B is a N x 1 matrix as well. In that formulation, the solution of a 
differential equation involves is the root of the function F = A ■ x — B. This 
means F(x new ) = 0, where new means the matrix x after an iteration. Thus 
F(x new ) = F(x id) + F'(x i d )(x new - x oU ) + ■ • • , which is in fact a Taylor 
expansion where terms beyond linear are considered unimportant. Solving 
in terms of 8 = x new — x Q id, we have: 

5 = _Fk%m) ( 4) 

F'(x old ) 1 ; 

where the prime denotes differentiation with respect to x. The above equa- 
tion is the correction to our initial guess x and iteratively leads us to the 
solution. Practically, one does the following steps: 

(i) Write the equation. In our case we have second-order differential 
equations. The derivatives are replaced in the following way: 

x , x (8.15) 

(ii) Write matrices A and B. In this case, matrix A is tri-diagonal. B is 
initially zero unless the equation is non- homogeneous with respect to x. 
Otherwise, it contains these non-homogeneous terms. If the equation 
is non-linear, then A = Al + A^l, where L and NL stands for Linear 
an Non-Linear respectively. 

(iii) Compute the inverse of the Jacobian of A. A general formula for 
the calculation of the Jacobian of a matrix is: 

dF 

J„ = ^ (8.16) 

The Jacobian, and its inversion can also be done at once with clapack 
1 routine package in order to achieve greater efficiency 2 . 



More information at www.netlib.org 

2 In fact it was observed that by using this package, the program needed at most 1/3 
of the time it needed before in order to run. The difference was greater by increasing the 
grid points. 
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(iv) Compute correction. The above step does the correction and this is 
added to our initial guess in the first iteration. Thus, our initial guess 
is improved in every iteration. 

(v) Find solution. When the correction S becomes smaller than a user- 
defined small value, the program stops and we get the numerical solu- 
tion. 

• Concluding remarks for Newton-Raphson 

This method, as one can observe in Part I, is one of the two methods which 
gives satisfactory results on both models described there. A result is sat- 
isfactory when we have proofs 3 that we found the solution, such as the 
appropriate asymptotics, the satisfaction of the equations of motion/field 
equations and of virial theorem. Apart from this, other advantages of this 
algorithm is the relatively short time it took us to get the results shown for 
both "test-models" and it is also the fact that it is relatively easy (compared 
with Relaxation and Runge-Kutta) to use. 

8.4 Energy minimization 

• How it works 

The general idea is this. Suppose the energy functional of a model of interest. 
This can involve two or more functions of one or more variables. One can start 
with an initial guess for these functions approximating the true solution. A 
standard minimization algorithm such as the one we used, adds a correction 
to the fields but having as a criterion that every such step will reduce the 
energy functional value. Thus, it is called energy minimization. The output 
is a field configuration for which the energy functional value is minimum. 
The algorithm needs as input the function as well as the derivative of that 
function at every point. 

Details of the specific algorithm can be found on page 428 of [104] . There 
are used two subroutines, the first of them called DFPMIN which has as 
target to find the so-called "Newton" direction which is a direction that 
makes the value of the energy smaller. Then, after this direction is calculated, 
subroutine LNSRCH is called and finds an appropriate step of correction 

3 In Part I, one can find in detail, explanation of ways of how to check if a final configu- 
ration is the solution of a system of differential equations, with the help of virial relations. 
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for every point of every function on the "Newton" direction. When the 
correction is done, DFPMIN calculates a new "Newton" direction based on 
the new corrected data. When the corrections are smaller than the machine 
tolerance (~ 10~ 8 ) or the derivative of our function is very small then we are 
on a minimum. 

To determine the next iteration point in this algorithm we do the follow- 
ing: Consider finding a minimum by using Newton's method to search for a 
zero on the gradient of the function. Near the current point Xj we have to 
second order: 

/(x) = /(xi) + (x - X() ■ V/(xi) + l(x - Xi) • A • (x - x,) (8.17) 

so 

V/(x) = V/(x i )+A-(x-x i ) (8.18) 

In Newton's method we set V/(x) = to determine the next iteration step 
and finally we have: 

x m = x, + A- 1 ■ (V/(xi +1 ) - V/(x,)) (8.19) 

The DFPMIN subroutine provides the appropriate formulas for calculating 
A -1 . One can find more details on p. 428 of [104] . 

• Numerical details 

A point we need to mention here is that one has to be very careful when 
inserting the derivatives (below denoted as gj) of functions (below denoted 
as pi) such as the one we have in our model. One must write down explicitly 
the first two or three terms of the double sum which constitutes our energy 
functional, as well as the last two in order to see the sequence of appearance 
of pi in order to compute the relevant g^. Then, we choose the variable of 
differentiation and gather the terms which include that variable in order to 
find the derivative. In numerical language, the notion of the term "variable" 
is not necessarily the same as in mathematical language. 

In this model, as we have five functions of two variables, we firstly have 
to divide the p-coordinate which has n points, into five parts representing the 
functions: 1 — > k\ for F, k\ + 1 — > k 2 for P, k 2 + 1 — > k 3 for A v , k 3 + 1 — > fc 4 
for W p and k± + 1 — > n for W z . We also have to divide the ^-coordinate in, 
say, n z points. Because of the fact that now we have derivatives of these 
functions also in the ^-direction we have to find a way to connect a point of a 
function in zu level with the same point in the zu + level. Well, according to 
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our formulation, this happens when we add n to the index of a point. Thus, 
if i is for p and ii is for z, then we can represent a derivative of a function / 
in z direction as 

df _ f[i + n]-f\i] 
dz 



with 



h. 



n 7 . 



(8.20) 



f8.2i: 



Also, differentiation of / with respect to p is 

df f[i+l]-f[i\ 



with 



dp 



h, 



Pmax Pr. 



n 



[8.22) 



(8.23) 



In the formulation below, for simplicity, we represent with ii + the next z 
level and with ii_ the previous z level. For example if we are at z = d level 
then ii — d-n and ii + = (d+ 1) • n, while u_ = {d — 1) • n, where d an integer. 
The energy functional in numerical language, follows: 



Ei:l^ki-l, ii: 0^n(z max -l) = 2l{V\h p h z pi 
Pk 2 +i+ii+ Pk 2 +i+ii 

h. 



1 

2n 



Pk2+i+l+ii ~~ Vk 



2 + '+" 



+ 



|_ 1 j Pk^+i+l+ii ~ Pk^+i+ii 
hn 



Pk 3 +i+ii + Pk 3 +i+ii 

h* 



+ 



Pki+i+l+ii Pki+i+ii 



Pkx+i+ii^ Pki+i+ii 

h* 



Pi+l+ii ~ Pi+ii 

hn 



+ 



Pi+ii + Pi+ii \ Pk±+i+ii 



h. 



+ 



Pi 



ePk 2 + 



N) + 



+ 



qp k3+ i +ii - M 



9i ( 2 



zacos 2 Q 

(Pi ~ a f 



qp ki +i+ii + M 



cos 2 6 
Pi - a 



Pi+ii 



+ a \Pi+i; 



92 f 2 



4 



+i+ii 



— U 



93 2 2 
r, Pi+iiPk!+i+ a 
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Below we write down the derivatives of the above energy functional for every 
"variable". We start with the i — 1 level for all functions and all ^-levels 
except for the first and the last one. In fact i = 1 means p = 0. 



Si+ii = 27w 1 h p h z p 1 
+ 2pi+u 



2 

— - Pl +ti _) - J^(P2+H - Pl+ii) - Jg[Pl+ii+ - Pl+ii) 



h 2 



qp ks+1+ii - M 



zucos 2 Q 

(pi - a Y 



qpk i+ i+u + M 



cos 2 Q 
Pi - a 



+ 



+ 9iipi + u - + 93Pi+iiPl 1+ i- 



Jki+i+H = 2irv 1 h p h z p 1 



TpiPki+l+ii ~ Pkx+l+ii-) ~ J^iPkx+2+ii ~ Pki+l+ii) ~ 

2 

_ J^\Pki+l+H+ — Pkx+l+ii) + 



+ ^Pkx+i+iiiepkz+i+u + N) 2 + g 2 (pl 1+1+u - u 2 )p kl+1+ii + g 3 p 2 1+ nPki+i- 
Pi 



gk 2 +i+ii = 2irv 1 h p h z p 1 



7^2 ( ~ J^(Pk2+2+ii — Pk 2 +l+ii) ~ 
2 

_ J^\Pk2 + l+H+ ~ Pk 2 +l+ii) + 



+ J~2 \Pk2+l+ii ~ Pk 2 +1- 



l+l+M 



k 2 +l+ii 



+ N) 



?k 3 +i+ii = 2-Kv l h p h z pi 



1 / Pk 4 +2+ii — Pki+l+ii Pk 3 +l+ii+ ~ Pk 3 +l+ii 



h, 



h z 

1 I Pki+2+ii- ~ Pki+l+ii- Pk :i +l+ii ~ Pk 3 +l+ii- 



h, 



h. 



+ ZqpI+u qpk 3 +i+u - m 
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The derivatives of all functions on all the p and z levels apart from the first 
and the last level follow. 



g(i:2-ki-l) i+ii = ZlTVihphzPi-! 



(Pi+ii ~ Pi-l+ii) 



2nvih p h z pi 



2 2 2 

jr;(Pi+ii — Pi+ii-) ~ J^iPi+l+H ~ Pi+ii) ~ J^\Pi+ii+ ~ Pi+ii) + 



qp k3 +i+ii - M 



zacos 2 Q 
(pi - a) 2 



qpk 4 +i+u + M 



cos 2 Q 
Pi - a 



2-r 



+ 



+ Qiipt+u ~ l)Pi+u + QsPi+uPh- 



■i+n 



g(i: 2->ki-l), ki+i+ii = 2nVih p h z Pi-i 

+ 2irv 1 h p h z p i 



J^iPkt+i+ii — Pki+i-l- 



2 2 

TpiPkr+i+ii ~ Pkt+i+iiJ) ~ J^(Pk!+i+l+ii — Pki+i+ii) ~ 

2 

_ ~L2\Pki+i+ii+ ~ Pki+i+ii) + 



+ ^P^+i+iiiep^+i+a + N) 2 + g 2 (pl 1+i+ii - u 2 ^p kl+i+ii + g 3 p 2 +ii p kl+i+ii 
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g(i: 2->ki-l), k 3 +i+ii = 1^V\h p h z pi 



1 | Pki+i+l+ii ~ Pki+i+ii 

h z y hp 
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2-KVihphzPi 
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cos 2 Q 
Pi - a 



The derivatives of all functions on the last p level and all z levels apart from 
the first and the last level follow. 



gk 1+ ii = 2nv l h p h z p kl _ l 

gk 2 +ii = 2nv 1 h p h z p kl _ 1 
1 



li 2 



[Pki+ii ~ Pki-l+ii) 



{Pk 2 +ii ~ Pk 2 -l+ii] 



gk 3 +ii = 2irv 1 hph z pk 1 - 



{Pk 3 +ii — Pk 3 -l+ii) 



gk 4 +ii = o 



gn+ii = 27TV l h z p kl ^i 



Pn+ii Pn—l+ii Pk&— Pfc4— 1+i 



h, 



h- 



135 



Summary of the algorithms used 



The derivatives of all functions on all the p levels apart from the first and the 
last level and the ii — level follow. The ii = level is in fact the z = —z . 



g(i: i+0 = 2-KVxhph z pi-l 



— (p i+0 -Pi-l+o) 



+ 



+ 2irvih p h z pi 



2 2 
j^(Pi+a + -Pi+o) - -^(Pi+i+o - Pi+o) 



Mz cos 2 e 
(Pi - a) 2 



)'*( 



M cos 2 e 
Pi - a 



9iPi+o (p 2 i+0 ~ l) + g3Pi+oPl 1+i+0 



+ 2nvih p h z pi 
2pk!+i+o 



2^ki-i), ki+i+o — 2nvihph z pi-i 
2 



(Pfci+i+O — Pfei+i-l+o) 



+ 



^ 2 (Pfci+i+1+0 ~~ Pk ± +i+o) — J^{Pk 1 +i+ii + — Pfci+j+o) + 



+ — ^ — ( e Pk2+i+o + N) + g 2 Pk 1+ i+o [Pk 1+ i+o ~ u ) + 93P l+0 Pk 1+ i+o 



g(i:2^k!-i), k 2 +i+o = 2nvih p h z p 



i-i 



2p 2 -i V K 



i.Pk 2 +i+0 — Pk 2 +i-l+o) 



+ 



2ixv l h p h z p l 



2pf K 



{Pk 2 +i+l+0 ~ Pk 2 +i+o) ~ J-^{Pk 2 +i+ii+ ~ Pk 2 +i+o) ) 



+ 



2<_ 



■j+0 



(ep k2+i+0 + N) 



g(i: 2-^-1), k 3 +i+o = 2-rrvxhphzp 

Pk 3 +i+ii + ~ Pk 3 +i+0 \ 



J_ / Pkj+i+l+0 — Pkj+i+0 
h z \ h p 



Mz cos 2 e 
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8.4. Energy minimization 



2->ki-l), k 4 +i+0 = ^Vih p h z pi-i 



1 / Pkj+i+O — Pkj+i-l+O 

h p y h p 

Pk 3 +i-l+ii + — Pk 3 +i-l+ 

h* 



+ 2irv 1 h p h z p l 



1 / Pfc 4 +»+l+Q ~ Pkj+i+O _ Pk 3 +i+ii+ — Pk 3 +i+0 \ 

hp y hp h z J 



+ 



+ 2 QP 2 i+o \qPk i+ i+o + 



M cos 2 e 
Pi — a 



The derivatives of all functions on all the p levels except for the last level 
and the last z level follow. The last z level ii max is in fact the z = zq. 



l->ki-l), i+iir, 



2nvih p h z pi 
2 



hi 



iPi+iimax Pi+iimax-l) 
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1 / 2 



hi 
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ma x 
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Summary of the algorithms used 
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Pmin Pmax 

F(p, z) 


Pmin Pmax 


Pmin Pmax 

A v (p,z) 


Pmin Pmax 

W p (p,z) 


Pmin Pmax 

W 2 (p,z) 



z = -2 



Above one can see how a grid, like the one we use, looks like. The five 
areas where the data points of every function are stored, is clear as well as the 
p and z levels. A brief diagram of the program we use follows. Subroutines 
DFPMIN and LNSRCH are provided by [104] . In the following diagram x 
is the matrix that has all the points of the five functions and is the matrix 
on which the corrections are applied until we have the final position of these 
points, the final configuration for the fields and the value of the total energy. 
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8.4. Energy minimization 



1. MAIN PROGRAM 

► Insert parameters, initial guess 

► Compute initial E and V-E" 

► CALL DFPMIN-Compute Newton direction p 

► Get final configuration from DFPMIN -< 

► Compute virial, Produce plots etc. 

► END PROGRAM 



2. DFPMIN 

► Insert data from main 

► Compute Newton direction p on initial guess 

► CALL LNSRCH-Compute step on p and correct x 

► Get corrected x from LNSRCH 

► IF Ax OR VE < 5 • KT 8 RETURN SOLUTION - 

► ELSE Compute new Newton direction p 

► END DFPMIN 



3. LNSRCH 

► Insert data from DFPMIN 

► Start with full Newton step A = 1 

► Correction: x neiu = x oM + Ap 

► Compute E new (x new ) 

► IF A < 5 ■ KT 8 OR E new < E old + Xag' RETURN 

► ELSE IF A = 1 THEN X new = -g /2{E new - E old - g) 

► ELSE IF b 2 - 3ag(0) < THEN X new = | , (a, b given) 

► ELSE X new = (-b + v^ 2 -3ag'(0))/3a 

► IF X new > 0.5A THEN X new = 0.5A 

► IF X new < 0.1 A THEN X new = 0.1 A 

^ A A neu! , E id E new 

► END LNSRCH 



4. FUNC 

► Compute E 

► RETURN the value of E 

► END FUNC 



5. DFUNC 

► Compute the derivatives of E 

► RETURN the values of the derivatives of E 

► END DFUNC 
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Summary of the algorithms used 



• Concluding remarks for Energy minimization algorithm 

It was an algorithm which successfully generated the expected results for our 
test models in the shortest time and the output for every model is exhibited 
in the appropriate section. It is the easiest to handle, even when compared 
to Newton-Raphson, and the results again were multi-tested as one can see 
in detail. Thus, it is used in our research models and provides the solutions 
we found there. 
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1.2 Vortex-antivortex formation in a nematic liquid crystal. Further 



1.3 Higgs field space and physical space. Above T c the potential has 
a minimum at <j)\ = = fa, but below T c the new form of the 
potential has its minima on a circle. If m is an element of the coset 
space G/H and \&i and ^2 two different states at that circle, then 
m^i = *$>2- Generally, every element of the coset space G/H can 
produce all the other minimum energy states when acting on any of 
them. Thus, G/H = M where M stands for the vacuum manifold. 
On the right side of the figure one can have a picture of how a 
cosmic string forms when different regions of space met, having 
their Higgs field phases as above. Picture from astro-ph/0303504 



1.5 Temperature dependent potential. Here is an example of such 
potential and how it changes shape as temperature falls. It is 
clearly visible that above T c the system has one stable minimum 
at <j) c = 0. As T — > the potential develops a minimum at (f) c 7^ 
and when the system settles there, we say that the initial symmetry 
of the system is broken. The ground state of the model does not 
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tails can be found in [TTJ 1 1 5 



details can be found in [T2] 16 



(A.Gangui) H7| 



1.4 Cosmic string formation. Picture from|astro-ph/0303504 (A.Gangui). 1 1 8 



astro-ph/0303504 (A.Gangui) Hi 



1.6 Apart from the intersection of two different strings (a), self-intersection 
of a cosmic string can be another mechanism of loop production 



(b). Picture from [135] [22 
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2.2 An example of a second order phase transition. The plots are for 
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2.4 Nielsen-Olesen string. The graph on top is the output of a relax- 
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Raphson algorithm and the bottom graph is done through energy 



minimization. All graphs show f(p) and a(p), for n = 1, A = 2e 2 = 4. |42 



2.5 Straight bosonic superconducting string. Both graphs on top are 
the output of energy minimization algorithm. Parameters are 
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a field 4 



3.1 Goldstone model: The solutions of equations (13.250 and (|3.26p 
are plotted as functions of k 2 for two values of (n, m). The dashed 



curve indicates the limit (|3.24[) [55 



3.2 Goldstone model: The energies of some of the solutions are plot- 
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branch indicate the four bifurcation points on the N = 2 type-I 
solution. The two lower stars indicate the two bifurcation points 
on the N = 1 type-I solution. The two numbers in parentheses re- 
fer to the number of negative and zero modes of the corresponding 



solution 58 



3.3 Goldstone model: The values (|3.1ip are plotted as functions of 
L for j = 0, 1, 2 (solid curves), together with the values cof(l, 1, — 1), 
wf (1, 0,-1), ^1(1, 0, 1), (4(1, 2,-1), lu 2 (1, 1, 1), cu 2 (l, 3, -1) of 4336]) 
(dashed and dotted curves). The numbers indicate the multiplicity 



of the eigenvalues [60 



4.1 U(1)a model: The relative position of the supercurrent as well 
as of the magnetic field on a xy-profile of the system on the right. 



The left picture is strongly reminiscent of an infinite solenoid. . . . [68 
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figEU) |73| 
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of P(p) and Ba{p) as a goes from 2.17 (dashed and dotted - • - ), 

to 2.15 (dashed - - - ) and 2.13 (dotted • • •) [74 
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5.1 U(1)a x U(l)w model: Step 1: Infinite straight string. Step 2: 
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tential (eq lfTTE]) for (gi,g 2 ,u) = (10,8, 1) [87 
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radius E(a). The middle graph is the supercurrent 1^ vs. a. In 
that graph one can clearly see current quenching. The bottom 
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radius E(a). The middle graph is the supercurrent 1^ vs. a. In 
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z = plane [94 



soliton we are after 94 



same in both of them thus, we choose to plot only the right side. . |95 



5.8 U(1)a x U(l)w model: Typical plot of the final configuration of 
fields. Parameters are M = 1, N = 1, e = 5, q = my/ = 2, 
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5.11 U(1)a x U(1)w model: The plot of the energy of the system 
vs. the radius of the torus for four different sets of parameters. 
Dotted is for (51,52,53)= (14,12.5,14), dashed is for (18,15,18), 
dashed and dotted is for (25, 20, 25) while dashed with two dots is 
for (30, 24, 30). For all data sets we have (e, M, N, u)=(10, 1, 1, 1). 
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5.12 U(l)w model: Profiles for the scalar field F(p,z) on the left and 
for the magnetic field By/ on the right side, for three random values 
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